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SUMMARY: When a flow with a non-uniform velocity passes through a cascade, 
the blades impose on each layer of the flow local pressure changes which are 
roughly proportional to the dynamic head in the layer. This results in local 
pressure differences which, within the blade passages, produce rotational 
secondary flows and which also cause overall spanwise accelerations and 
displacements. These effects are analysed for frictionless flow by treating the 
blade row as an actuator plane whose outlet angle is determined by the 
secondary velocities in the blade passages. 


At low outlet angles or large deflections the secondary flows effectively 
control the downstream flow angles. When the outlet angle is large the 
displacement effect can result in downstream flow angles of opposite sign to 
those induced by the secondary flow alone. 


Measurements on cascades with small spanwise variations induced in the 
approaching stream confirm the main predictions of the first order theory. 


31. Introduction 


Flow with a non-uniform velocity through a cascade of aerofoils is of interest 
in the theory of axial-flow turbo-machines because the non-uniform flow in the 
boundary layer region near the annular walls causes a loss in efficiency and a 
reduction of the effective area for flow. 


In some theoretical”) approaches to the study of three-dimensional flow 
in turbo-machinery use has been made of the actuator disc concept in which the 
blade row is replaced by an infinitely thin axially symmetrical disc which, as in 
propeller theory, imposes a sudden change. More detailed studies®’™® of the 
behaviour of the fluid in the passages between the blades have shown that, in 
regions of initially non-uniform flow, secondary velocities normal to the main 
direction of flow are generated, with consequent changes in the lift and circulation 
about the blades. 


This paper describes a theoretical and experimental study of the flow through 
a cascade of untwisted aerofoils placed in a wind tunnel of rectangular cross 
section in which the approach velocity was varied artificially in the spanwise 
direction only. The variation in velocity was obtained by placing an obstruction 
upstream, which gave a two-dimensional wake which was symmetrical about the 
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centre of the blade span. The velocity in the wake varied approximately 
sinusoidally in the spanwise direction. 


The non-uniform flow in the wake produced secondary flows in the blade 
passages which altered the circulation about the blades, and caused the angle of 
flow leaving the blades to depart from the angle usually measured in uniform flow. 
As the secondary flows were not large, the theory of Ref. 4 could be used to 
predict the variation in outlet angle. 


The row of blades was then imagined to be an actuator plane, which produced 
a deflection of the flow and which had an outlet angle at the trailing edge position 
given by the two-dimensional value modified to take into account the effect of 
the secondary flow in the blade passages. This enabled the flow angles and 
velocities far downstream to be determined and compared with the results of the 
experiments. 


The theory is developed for inviscid, incompressible flow and the experiments 
were performed with turbine or accelerating cascades to reduce the possibility of 
boundary layer separation on the blades. 


NOTATION 
a_ velocity decrement defined by equation (23) 


C constant 
K=1-exp(- =x./8) 
m=t?/(2+?’) 
m=(14+K)/[2+(1+ 
p, stagnation pressure 
s blade spacing 
=SCOS 
t= tan 2,, (independent of z) 
t/=t—tan a,, (a function of z) 
t”=t- tana, 
u,v,w velocity components (see Fig. 1) 
u, value of u, outside the non-uniform region 
u, Value of u, outside the non-uniform region 
u, axial velocity at the trailing edge 
Ua=U—U, 
U,=u, seca, 
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f= CONSTANT 
x=0 x 


ACTUATOR PLANE AT TRAILING EDGE 


Fig. 1. 
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ACTUATOR PLANE NOT AT TRAILING EDGE 


Diagrams showing nomenclature. 


V_ velocity vector with components u, v, w 


V_ average value of velocity V in the Y-direction 


X,Y,Z 


Xe 


Cartesian co-ordinates (see Fig. 1) 


distance from actuator plane location to the trailing edge 


Y co-ordinate normal to the undisturbed stream at outlet 


Goo 


€ 


¥ 
X 


$= 


upstream angle of approach to the cascade 


the angle at which the flow leaves the trailing edge 


the angle at which the flow leaves the trailing edge outside the non- 


uniform region 

4 (boundary layer thickness) 
deflection (2,, — ,) 

density 


stream function 


stream function of the secondary motion 


curl V, vorticity vector with components ¢ 


Suffixes 


0 
1 


vorticity components 


refers to conditions at the actuator plane 


refers to conditions upstream of the cascade 


2. refers to conditions downstream of the cascade 
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2. Actuator Plane Theory 


2.1. GENERAL SOLUTION 


The nomenclature and co-ordinate system used is shown in Fig. 1. The flow 
is separated into two regimes on either side of the actuator plane. In each regime 
the equations of motion and continuity may be written in the form 


and div V=0 ‘ (2) 


where p, is the stagnation pressure, p the density, V the velocity vector with 
components u, v and w and €2=curl V, the vorticity vector, with components 
£, », ¢. These velocity and vorticity components correspond to the Cartesian 
co-ordinates x, y, z, where z is the spanwise direction (see Fig. 1). If the cascade 
is assumed to be many-bladed and therefore infinite in the y— direction, variations 
in this direction will not occur and the components of vorticity may be written 


The continuity equation becomes 


ou (5) 


It is therefore possible to select a stream function ¥(x, z) which is independent 
of y and which enables the intercepts of the stream surface, Y/=constant, on any 
plane y=constant to be determined. From equation (5), 


(6) 
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' Since the pressure is independent of y, the second Euler equation becomes 


ov ov 


whence v is constant along a stream surface intercept, Y=constant. 


The second Helmholtz equation yields 


On On dv 
Ha 


Ox oz (8) 


and, on substitution for € and ¢ from equation (3), the right hand side equates to 


' zero. Hence » is also constant for constant y. 


The solution of the equations of motion for each regime is generally 


_ intractable without some approximation. If the non-uniformity of velocity in the 


approaching stream is also small, then terms such as du/dz, dv/dz 
and » are small, and therefore the spanwise velocity or w/u is also small. Since, 
from equations (3) and (7), 


ov w dv 


{ and the variation of v with x are the product of two small terms and may be 
neglected. 


Similarly from equation (8) 


and the variation of » with x may also be neglected. This type of approximation 
was suggested by Ruden® for annular cascades. An exact analysis) of the flow 
through an annular cascade confirms that it is reasonable. 


Since the vorticity components are now, to the first approximation, independent 
of x, physical consideration of the axial velocity field induced by the »— components 
of vorticity on both sides of the actuator plane leads to the result that at the 
plane, x=0, the axial velocity is the mean of the values at x= — 00 and +00. 


Across the actuator plane, by continuity, 


Since the flow is inviscid, p, is unchanged, and since there are no blade forces 
in the z— direction, 
Wo1 =Wo2=Wo- ‘ ‘ & 
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The second of these conditions and the assumptions for » and v may be applied 
to equation (4) to give 


dv, ov 


where the values at x= +00 may be used for 7, and 7, and v, is obtained from the 
known value of the upstream angle of approach to the cascade, 2,, i.e. v, =u, tan 2, 
To obtain v,, knowledge of the outlet angle at the trailing elge of the blades is 
assumed. 


There is some evidence to show that the actuator plane should be located at 
the centre of pressure of the blading. The equations are somewhat simpler when 
it is located at the trailing edge and this position will be assumed first. The 
modifications required in the analysis for other locations are given in the Appendix. 
With the actuator plane at the trailing edge, v, =u, tan 2,,, where @,, is the angle 
at which the flow leaves the trailing edge. 


Substituting in equation (12), 


Uy +u, tana, (u, tan z,)=u, +u, tan (u,tana,.) . (13) 


Taken with equation (9), this gives a non-linear first order differential equation for 
u, which may be solved numerically or graphically. 


The equation shows that, if z, and 2,, are independent of z, then outside the 
region in which u, varies with z, uv, and u, are constant. This implies that the width 
of the non-uniform region, i.e. thickness of wake or boundary layer, does not 
increase, and is a consequence of the approximations made in deriving equation (12). 


Equation (13) may be linearised by writing u,—u,=2u’ and assuming that 
u’ is small compared to u,, so that terms of order uw’? may be neglected. Then 


du’ 
u, (2+ tan? +u de (u, tan? z,,.)=u, tan 2, (u, tan 


— u, tan 2,, (u,tana,.). . (14) 


With tan z, and tan2,, constant this may be solved to give 


(15) 


where m=??/(2+1?) and t=tan2,,, and C is a constant which may be found by 
applying the continuity condition. 
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Now v,=u, tan 2),=u, tan @,, and with the same degree of approximation 


tan a,,—tan ES COS? a5 
tan COS? 


1 | « @ 


2.2. APPROXIMATE SOLUTION FOR 2, CONSTANT AND %. VARYING SLIGHTLY 


When «g,, varies, an approximate solution may be obtained by writing 
tan 2,.=t+¢’, where ¢ is independent of z and ? varies with z. Assuming that ¢ is 
small, so that terms in ¢’? and u’t’ may be neglected, equation (14) becomes, with 
tan 2 constant, 


ut du, _ 2du,_ t d(u,*t’) 


If tan 2, may be written tan z,=f+?”, 
then, since V,=U, tan tan @,, 
(u,+u’) (t+ +2u’) (t+?t’). 
Whence, neglecting terms in u’f and u’t’, 
ut=u, (t—t” 
Substituting for u’ in equation (17), 


du, m du,  2u,t du,” d(u,t’) 


d. 
dz (u,t”)+ mt 


If outside the non-uniform region u, reaches a constant value uv, and ¢ —> 0, the 
solution is 


(20) 


omitting an arbitrary constant. 


The last term, which is in general small, may be further simplified by noting 
that, if the variation in u, is small, 
=[1-m(1- 
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Hence the equation becomes 


%,/ 

omitting an arbitrary constant. , 


The expression for u’ in equation (17) may be similarly treated to give 


2mu,t 2m (1 re ‘a (1 


where u, is the value of u, far outside the non-uniform region where it is assumed 
that *¢—>0. The value of u, is determined from continuity. Outside the region 
of varying u,, u.-u,=2miu, t/t, so that u, will rise or fall to its final value beyond 
the non-uniform region, depending on the sign of /. 


2.3. THE EFFECT OF SECONDARY FLOW IN THE BLADES 


The secondary flow in the passages between the blades alters the value of <,, 
from that found in two-dimensional flow. The variation in ,, is calculated by the 
method suggested in Ref. 4. The calculations are performed for a non-uniform 
flow which varies sinusoidally within a “boundary layer” region. The velocity 
within the boundary layer region of thickness 28 is given by 


1 1 +008 . ; ~ 


The stream function x of the secondary motion within the blade passages is given by 


2y 

where ¢ is the deflection (a,,—,), U, the approach velocity (=u, sec 2,), and Y is 
the co-ordinate normal to the unperturbed direction of the stream at outlet, i.e. a 
direction making an angle z,, with the x-axis (Fig. 1). The boundary conditions 
are x=0 at Y=0 and s’, where s’=scos @,, and s is the blade spacing, and x=0 at 
z=0 and z=00. The latter boundary condition implies that the region of uniform 
flow outside the “ boundary layer ” region is large. 


a solution is of the form X= Xn (z) sin nt Ake . (26) 


254 The Aeronautical Quarterly 


This 
uy ‘) uy 
and 
Whe 
The 
and 
The 


SHEAR FLOW THROUGH A CASCADE 


| This satisfies the boundary conditions at Y =0 and s’ and gives, 


and when z > 8, Xn (5) Xn=0. . (27b) 


When solutions to these two equations for x, and x,’ are matched at z=6:— 


for 0<z<, 


0 
and for z> 4, 


3 
(nny? e aA sinh (2 ) aa. . (28D) 

0 


| These equations satisfy the boundary conditions at z=0 and 00. Now the velocity 
V in the Y —direction is 0x/0z. Hence the average value of V at any z is given by 


(29) 


and the maximum value by 


n+1 


> (- 1)? 


and the corresponding change in average outlet angle (Fig. 1) by 
V 2 COS Xn’ 


Aa 


02 


The change in maximum outlet angle 


cos nth 
AQ — > (- 1) Xn 


Now = A tan 2,, ~ sec?z,,Aa,, 


2 SCC Bos Xn 
Tu, n 
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Using the value of u, from equation (23) and differentiating equation (28a): — 


For 0<z< 4, 


Hence, for 0<z< 


1 


nd 
RZ 
u a 5 
r(“) = 16¢ sec %,, sec x, (*) (33a) 
u 2 
1 8, 2a 
L Ss’ 


and for z> 4, 


a 
t/= — 16¢ sec 2,, seca, (2) 


In general, terms above n=1 in these series are small enough to be neglected. 


2.4. APPROXIMATE SOLUTION FOR DOWNSTREAM ANGLE 


Equations (33a) and (33d), together with equation (21), enable the variation of 
angle far downstream to be determined. The general behaviour is first shown by 
an approximate analysis. The integral term in equation (21) is neglected, as are 
terms above m=1 in equation (33a). 


If a/u, is small, then from equation (23), 


34) 
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Therefore, for z < 4, 


cos? Goo a TZ 
cos? a, 1) (m+1) + 16 (m+ 1) G2 - a,)x 


TZ 
a a COS +e 
X SCC SEC a, (£)(14 1+c0s =) 
uy, u, 


2.5. COMPRESSOR CASCADE 


For a compressor cascade | @,, |< |, |, and for a turbine | z,.| >|, |, so 
that the signs of both terms on the right hand side of equation (35) are reversed as 
the cascade changes from a compressor to a turbine. For a compressor cascade the 
first term, representing the effect of the spanwise displacement which occurs through- 
out the field of flow, is positive for all values of z in the non-uniform region. The 
second term is negative for small values of z, is positive at the boundary where 
z=6 and passes through zero at the point where cos (™z/5)+e~**'* cosh (z/s’)=0. 
It is possible therefore, to have either positive or negative values of ¢”, i.e. the 
change in tan 2, from some constant value well outside the non-uniform region, at 
and near z=0, but the value will always be positive at z=6. A positive value of 
t’, i.e. an increase of 2, above the main stream value, implies an “ under-turning ” 
or decrease in deflection, a negative value an increase. At z=0 the condition that 
t’ is positive and hence that there are no negative values of ¢” and there is no 


increase in deflection at any z is, from equation (35), 


SiN COS [ 1| 
16 (a,—%,) Lcos?a, 


2.6. TURBINE CASCADE 


For a turbine cascade, 
signs of ¢” and the terms in equation (35) are concerned. However, a positive value 
of ¢” implies “ over-turning ” or an increase in deflection, and expression (36) repre- 
sents the condition that t” is always negative and there is no increase in deflection 
at any z. 


For both turbine and compressor cascades the secondary flow tends to increase 
the deflection at small values of z. 


2.7. NUMERICAL ILLUSTRATIONS 


To illustrate the effect on angle, a number of calculations have been made, 
based on the sinusoidal velocity distribution of equation (23) and on equations (21) 
and (33). Fig. 2 shows results for a compressor cascade with z,=50°, 2,.=30° and 
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0-6 10 12 


4X, o,CHANGE IN ANGLE AT TRAILING EDGE 


DUE TO SECONDARY FLOW 


0-6 +o % 12 


At, CHANGE IN ANGLE FAR DOWNSTREAM 
—-—--INTEGRAL TERM NEGLECTED 


INCREASED 
DEFLECTION 


Fig. 2. 
Downstream angles for compressor cascade. 
a/u,=0°25, a, =50° =30°, 5/s’=0°5, 8/s=0°433. 


3/s’=0-5, for a velocity profile for which the value of u, at z=0 is half the main 
stream value, i.e. a/u,=}. The figure shows first the effect due to the spanwise 
displacement only (i.e. with  =0). There is an increase of angle, i.e. a decrease of 
deflection in the wake. Next is shown the variation in angle at the trailing edge 
due to the secondary flow in the blade passage. There is first an increase in 
deflection at the centre of the wake which is reduced until a decrease in deflection is 
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Ao<,,DISPLACEMENT EFFECT, t':0 


Act, CHANGE IN ANGLE AT TRAILING EDGE 
DUE TO SECONDARY FLOW 


Ax,, CHANGE IN ANGLE FAR DOWNSTREAM 
---—-— INTEGRAL TERM NEGLECTED 


to 


- 
ra) 
2 


Fig. 3. 


Downstream angles for two turbine cascades of 20° deflection. 
a, =30°, a/u, =0°25, 8/s’=0°5, 8/s=0-433 
a, =40° =60°, a/u,=0-25, 8/s’=0°5, 8/s=0-25 


’ Zoe 


» Zoo 


noted at the edge of the wake. The combined effects which are not purely additive 
(see equation (21)) are finally shown in a curve at the bottom of Fig. 2. The slight 
effect of omitting the integral term in equation (21) is also shown. 


Figure 3 shows the results for a very low deflection turbine cascade. Two sets 
of curves are shown; one for an outlet angle of 30°, for which the change of angle 
due to the secondary flow appears to dominate the variation of 2, far downstream; 
in the other, for z,.—60°, the displacement effects become more important. In 
both cases the deflection is 20°. In Figs. 2 and 3 the effect of neglecting the integral 
term in equation (21) is seen to be slight, and it is omitted in subsequent figures. 
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0-8 


COMPRESSOR 


‘ /o,= 0-25 Effect of velocity profile. 


a/o,+0-125 


Ae<., CHANGE OF ANGLE FAR DOWNSTREAM 
(NTeGRAL TERM NEGLECTED 


TURBINE COMPRESSOR 
30 
Fig. 5. 
Effect of boundary layer thickness 
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Fig. 6. 
Effect of outlet angle: compressor blades. - 


17 50°30, 3, =0:577 


220: 253-0-5 


Ax<,, CHANGE OF ANGLE FAR DOWNSTREAM 


Fig. 4 shows the effect of the magnitude of the velocity reduction in the wake. As 
might be expected from equation (35), the smaller is a/u,, i.e. the smaller the 
velocity variation in the wake, the less is the magnitude of the angle variations. 


The effect of boundary layer thickness is shown in Fig. 5, in which curves for 
8/s’'=4 and 1 are presented. As the 3/s’ term enters only into the expression for 7’, 
change in 5/s’ does not affect the spanwise displacement (the first term in equation 
(21)). The larger is 5/s’, the smaller is the value of the angle change ¢ due to 
secondary flow. For the turbine with z,,=60°, the reduction of the secondary flow 
eect between 5/s’=4 and 8/s’=1 is large enough to change the sign of Aca, at 
small values of z/5. The change of boundary layer thickness also has a large effect 
on the value of Aa, for the compressor cascade shown. 


Figure 6 shows the effect of varying the outlet angle of a compressor cascade 
at constant deflection. The conditions are similar to those in Fig. 3, except that 
the ratio of boundary layer thickness 5 to blade spacing s is constant (5/s=4), so 
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Fig. 7. 
Effect of deflection: turbine blades. 


~ € 
10° 30° 
0-433 
30 40 2.0.25 0-5 


-------a0" 6f 106) 


Ax: CHANGE OF ANGLE FAR DOWNSTREAM 


that as 2,, is varied from 10° to 50°, 8/s’ varies from 0-51 to 0-78 approximately. 
The increasingly large displacement effect with higher outlet angle is shown by the 
change in sign of Az, between the curves for z,.=30° and 2,,=50°. Comparison 
of Figs. 3 and 7 shows the large effect of deflection on the change in angle due to 
secondary flow. In Fig. 7 the large values for z,.=60° and «=100° are clearly 
beyond the limits of the analysis. 


In the curves in Figs. 2-7 the actuator plane is located at the trailing edge 
of the blades. This has been done to simplify the analysis. The modification to 
equation (21) necessary to allow approximately for other locations is given in the 
Appendix and curves showing the approximate effect of locating the actuator plane 
at the mid-chord position are given in Figs. 13-15. Moving the actuator plane 
forward from the trailing edge reduces the displacement effect, but does not alter 
the general inferences which have been drawn from the numerical results shown 
in Figs. 2-7. 
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/2.8. DOWNSTREAM VELOCITY PROFILE 


The axial velocity u, far downstream may either be obtained from equation (22) 


‘or from equation (18). In some cases when the span is large compared to the 


thickness 8 of the non-uniform region, u,=u,, which gives 
(1-2-9) 
With the approximations made in deriving equation (35), 


COS? Zoo ) a ( =) 
= (m+ 1+cos ) +, 


32 (@,.—2,) SiN a TZ [= (xz/8)+e-"’* cosh 


The resulting variation of u, is the same as that for constant 2,., equation (15) 
modified by an amount which is approximately proportional to 2mt’/t, i.e. the third 


‘term on the right hand side of equation (22). Since 2m/t=2t/(2+1*) this term will 


be zero for @,.=0 and 90°, and reaches a maximum of 1//2 at tana,=/2. So 


| that, over a range of values of t, u, will not be a strong function of ¢. At the edge 
| of the boundary layer z=8, t”—? will have the same sign as ?’, i.e. positive for a 
' compressor and negative for a turbine. Hence, at z=, u, will be less than u, for a 
' compressor cascade and greater for a turbine. 


At the centre of the wake, z=0, if the value of ¢ is small, tr” — ¢’ will be positive 


_ for a compressor and negative for a turbine. Hence, at z=0, u, will tend to be 
| less than uw, for a compressor and greater for a turbine. This conclusion remains 
approximately true provided that 


12, (COS? COS 2, 
g se cos?a, (a, = Zoo) sin (3/s’)"] (39) 


If this does not hold, u, will tend to be greater than u, at z=0 for a compressor 


and less for a turbine. 


Velocity profiles based on the assumption that u,=u, are shown in Figs. 2 
and 3. These profiles follow the general conclusions already given. As shown in 
the Appendix and in Figs. 13-15, the results are somewhat modified by any 
movement of the actuator disc from the trailing edge. 


i 3. Experimental Work 


| 3.1. ExperIMENTAL CONDITIONS 


The results of two sets of experiments from somewhat similar cascades are 


available for comparison with the theory. In both experiments a stream with an 


artificially produced velocity profile symmetrical about the mid-span position was 
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allowed to pass through the cascade, and measurements of velocity and angle wer 
made at distances far enough downstream to be equivalent to the downstream 
infinity location.* The first set of experiments described were made by W. D. 
Armstrong and the second by E. A. Pinsley’’. 


In each set of experiments, only small deflections were used and the blades 
were arranged as a cascade of turbine blades, so that there was a pressure drop 
through the cascade. In this way it was hoped to minimise losses. Three widely 
different outlet angles were used in each set of experiments and measurements of 
the cascade performance with a uniform inlet velocity were made in each case to 
establish the two-dimensional properties of the cascade. 


3.2. ARMSTRONG’S EXPERIMENTS 


In these tests,t a cascade was set up in the working section (5 in. x. 10 in.) of 
a pressure tunnel, supplied with air from a 14 h.p. fan. The major dimensions 
were: —blade chord 0-90 in., pitch 1-08 in. and blade span 10-00 in. The blades, 
of which seven to nine were used, were constructed on a 33° circular arc camber 
line using the profile shape given in Fig. 8. At the ends of the blades the cascade 
walls were continued downstream for several chords. Boundary layer control 
was applied to the tunnel walls lying parallel to the blade span. This enabled the 
main streamlines in the neighbourhood of the wall to be slightly re-orientated and 
thereby ensured a uniform flow onto all the blades. 


The upstream velocity profile was artificially modified by producing a “‘ wake” 
upstream of the mid-span position, this wake being shed from a perforated plate 
12 in. long which divided the tunnel into two passages, each 5 in. square. The 
plate was at least 15 in. upstream of the centre point of the cascade. 


In all experiments, measurements were made of the inlet angle 2,, which was 
uniform, and the velocity one chord upstream and the velocity and angles one 
chord downstream. The measured angles and velocity profiles are given in 
Figs. 10-12. The outlet angles presented are those measured one chord downstream 
at a point midway between two blades, i.e. the maximum and not the mean. 


3.3. PINSLEY’S EXPERIMENTS 


A similar set of experiments was done earliert by Pinsley’”. The cascade 
tunnel and experimental techniques were very similar, although in this case the 
blades had a lower camber and there were no walls downstream of the cascade. 
The dimensions were:—chord 0-80 in., pitch 0-80 in., parabolic camber 22°6° 
and blade span 10-00 in. Fig. 9 shows the blade arrangements used, while Figs. 
13-15 give details of the incident velocity profiles, the outlet velocity profiles and 


*As may be seen from examining equations (A5) and (A6) in the Appendix and the details 
of the configurations tested, this condition is practically reached at a distance of one chord 
downstream. 

tIn the Engineering Laboratory, Cambridge University. 


tIn the Gas Turbine Laboratory, Massachusetts Institute of Technology, U.S.A. 
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the values of the velocity V in the Y—direction. The last were taken 1} chords 
downstream from the trailing edges, but the axial velocity measurements were made 

| 4% chords downstream. The incident profiles were measured about 8 in. upstream 
of the cascade. 


scade 
e the 
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22:6" 
Figs. 
>and | 4 Comparison of Theory and Experiment 
— § 4.1. ARMSTRONG’s EXPERIMENTS 
details 
chord & Figure 10 shows the results of tests with a cascade with a nominal outlet angle 
| of 25-5°. The upper curve shows the inlet velocity profile U,. A cosine curve 
with the relation of equation (23) was fitted. In the calculation of the 
variation of velocity and angle downstream the actuator plane was located half a 
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Fig. 10 
Cambridge experiments. 
a,=11:2°, Reynolds 
number=2°72 x 104. 
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° EXPERIMENTAL POINTS 


blade width upstream of the trailing edge. No pressure plots were made on the 
blades, so that the exact location of the centre of pressure is not known. This 
appeared, however, to be a reasonable position and the theoretical results with the 
plane located at the trailing edge give too much emphasis to the displacement effect. 
Consequently equations (A 12) and (A 14) in the Appendix were used to determine 
the outlet angle and downstream velocity profile. Equation (A12) had to be 
slightly modified since, due to the difficulty of interpreting angle measurements in 
the wakes behind the blades, the angles measured were the maximum rather than 
the mean. The value of ¢’ used in equation (A 12) was therefore determined from 
equation (30a) and not equation (30), but in calculating u, from equation (A 14) 
the mean value of ¢ was substituted. 


The comparison between theory and experiment shown in Fig. 10 is good, and 
Fig. 11, for the cascade with a nominal outlet angle of 40°, shows similar 
satisfactory agreement. In Fig. 12, with an outlet angle of 54°, the agreement is 
still possibly satisfactory but better agreement could probably be obtained by 
placing the actuator plane slightly farther forward. Comparison of Figs. 10, 11 and 
12 confirms the general prediction of the theory and shows that quite different 
variations of outlet angle z, and apparent deflection (z,—«,) are obtained in the 
three configurations. 
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4.2. PINSLEY’S EXPERIMENTS 


Figure 13 shows the results of tests with a nominal outlet angle of 30° and 
Figs. 14 and 15 the results with outlet angles of 45° and 60° respectively. The 
inlet velocity profiles used in these experiments gave larger velocity gradients and 
therefore inlet vorticity. The deflections were smaller, being of the order of 8°, 
68° and 5° respectively compared with the values of 14-4°, 11° and 14-6° used 
in Armstrong’s experiments. In view of this difference and its likely effect on the 
centre of pressure, theoretical calculations have been made with the actuator plane 
at the mid-chord position (full lines) and at the trailing edge (dashed lines). The 
value of 5 has been rounded off at 2-0 in. and, while it is clear that this gives 
theoretical inlet velocity profiles slightly larger than the measured ones, the 
measuring station was so far upstream (about 8 in.) that some increase in thickness 
of this upstream wake is to be expected between the measuring point and the leading 
edge. The maximum values of V, the velocity in a plane of constant z in the Y - 
direction (Fig. 1), i.e. normal to the nominal direction of flow, have been compared 
with the values reported by Pinsley. To calculate V the maximum value of Az, 
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Fig. 12. 
Cambridge experiments. a, =39-6°, Reynolds number =4-15 x 


was used, together with Pinsley’s measured velocity. These values of V differ from 
Aa, only by a multiplying factor which is approximately constant. The agreement 
between theory and experiment shown in Fig. 13 is not as good as that shown in 
Fig. 10, although the character of the variation of outlet angle is very similar. 
The outlet angle is larger at the mid-span position than the value given by the 
full-line curve which is calculated with the actuator plane at the mid-chord position. 
Presumably if the actuator plane were placed at the leading edge better agreement 
would be obtained. The minimum of the experimental values occurs at z=3in., 
whereas for both theoretical curves it is at about z=3} in. The experimental 
downstream velocity profile is also shallower than the theoretical one. This may 
be due partly to the effect of friction and mixing in the 10-5 in. between upstream 
and downstream measuring positions, which is ignored in the theoretical treatment. 
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Fig. 13. 
M.LT. experiments. «,=22°, «,,=30°. 


5 Similar differences between theory and experiment are noted in Fig. 14, 
2, ~~ 45°, although the agreement for the value of V,,.x at the mid-span position is 
good when the actuator plane is placed at the mid-chord position. 


The results shown in Fig. 15, 2,, ~ 60°, show a better correlation between 
theory and experiment with the actuator plane at the trailing edge. This is possibly 
justified, since the incidence decreases from —10-9° at Z,-— 30° to —13-9° at 
60°. 
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5. Conclusions 


The theory predicts a variation of types of distribution of angle downstream 
which are, for instance, markedly different for differing cascade staggers. These 
predictions are largely confirmed by the experimental results. The correlation 


between theory and experiment is best when the vorticity initially present in the , 


approaching flow is small, as in Armstrong’s experiments. Reasonable correlation 
is achieved for Pinsley’s results. The position chosen for the actuator plane is 
important and the correlation between theory and experiment reflects any 
uncertainty in the location of this position. 


At small nominal outlet angles the secondary flow in the blade passages 
largely determines the variation in angle downstream. At large outlet angles the 
secondary flow effects are masked by the large spanwise displacements of the flow; 
this is accompanied by a marked smoothing of the axial velocity profile. 


Other effects being equal, a thicker region of non-uniform flow reduces the 
initial vorticity and therefore the departure of downstream flow angles from the 
nominal values. The same applies to the effect of reducing the initial velocity 
variations. The larger the deflection is, the larger is the secondary flow in the blade 
passages and its influence on the downstream flow. 


Although the analysis ignores losses and is only valid for initially small non- 
uniformities in velocity, it should be of some assistance in the interpretation of the 
behaviour of the boundary layers on the walls of cascades or compressors. 


ACKNOWLEDGMENTS 


Acknowledgment is due to the Chief Scientist, Ministry of Supply, for 
permission to publish the results obtained from experiments in the Engineering 
Department, Cambridge University. 


REFERENCES 


1. Marsie, F. E. The Flow of a Perfect Fluid through an Axial Turbomachine with 
Prescribed Blade Loading. Journal of the Aeronautical Sciences, Vol. 15, p. 473, 
August 1948. 


2. Braca, S. L. and HAWTHORNE, W. R. Some Exact Solutions of the Flow through Annular 
Cascade Actuator Discs. Journal of the Aeronautical Sciences, Vol. 17, p. 243, April 1950. 


3. Squire, H. B. and Winter, K. G. The Secondary Flow in a Cascade of Aerofoils in a 
Non-Uniform Stream. Journal of the Aeronautical Sciences, Vol. 18, p. 271, April 1951. 


4. HAwTHorNE, W. R. and ArMsTRONG, W. D. Rotational Flow Through Cascades. Part II. 
The Circulation about the Cascade. Quarterly Journal of Mechanics and Applied 
Mathematics. Vol. VIII, Part 3, p. 280, September 1955. 


5. RUDEN, F. Investigation of Single Stage Axial Flow Fans. N.A.C.A. T.M. 1062, 1944. 


6. ARMSTRONG, W. D. The Non-Uniform Flow of Air through Cascades of Blades. 
Cambridge University Ph.D. Thesis, 1954. 


7. PinsLey, E. A. An Investigation of Secondary Cascade Flows. Massachusetts Institute of 
Technology M.Sc. Thesis, 1950. 


272 


The Aeronautical Quarterly 


t 


stream 
These 
elation 
in the 
elation 
ane is 
S any 


ssages 
les the 
flow; 


es the 
m the 
locity 
blade 


non- 
of the 


y, for 
eering 


Blades. 
‘ute of 


luarterly 


SHEAR FLOW THROUGH A CASCADE 


Appendix 


EFFECT OF ACTUATOR PLANE LOCATION 


When the actuator plane is not located at the trailing edge it becomes necessary 
to estimate the value of the axial velocity u, at the trailing edge, since the theory 
assumes that the angle a,, is the value at the trailing edge and hence v,—v, tan a,,. 
From equation (3), the vorticity downstream of the actuator plane is given by 


Ou Ow Ou, Ow 
"= Ox oz Oz Ox’ 


where u=U, + Ug. 


Since » does not vary with x within the limits of the assumption made in the main 
text, 7 = 0u,/0z and hence 


0z 


~ Ox 
Ou, ow 


ox 


Also, by continuity, 


Hence 
aw 


and 


The boundary conditions are that w=0 at z=0 and on any other wall. However, 
taking into account that, as mentioned in the discussion after equation (13), the 
assumptions for » imply that u, and u, are constant outside the region in which u, 
varies with z, it seems reasonable to assume that at the edge of the boundary layer 
du,/0z=0. With these conditions 


where m= 1, 2, 3 and so on. 


Neglecting all terms in the series except that for m=1, the velocity at the trailing 
edge becomes 


U,=U, — (u, —u,) exp (- (A6) 


where x, is the distance from the actuator plane location to the trailing edge. 
Hence v,=u, tan (u,+u’+ Ku’, 


where K=1- exp (—7x,/8). 
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Neglecting terms in ¢’?, u’? and u’?’, 


dv 


v,— Zz + +2u,fu 4K t?u’). 


The additional term in K needs to be included in equation (17), giving 


au, t d 
2 2)—— 1 tan — — — —(y,7 


Similar rearrangement gives, for equation (19), 


d du, m (1- du x) 


(A 10) 


where m’ =?? (14+ K)/[2+(1+K) (All) 


Hence the solution for ¢” similar to equation (20) and omitting the integral term which 
has been shown to be small is 


1+(1+K) ?? cos? a, 24+(01+K)2 
(A 12) 
Putting u, tan ,,=u, tan a,, we obtain an equation analogous to equation (18) 


From which, assuming 
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The Design of Intermediate Vertical Stiffeners 


on Web Plates Subjected to Shear 


K. C. ROCKEY, Ph.D., M.Sc.(Eng.), A.M.I.C.E. 


(Engineering Department, University College of Swansea) 


SuMMARY: The correct design of intermediate vertical stiffeners on web plates 
subjected to shear becomes very important when the web plates are designed 
to operate at loads close to their buckling loads. This paper presents details 
of an extensive series of tests conducted on stiffened web plates subjected to 
shear. From the analysis of the results obtained from these tests, new empirical 
relationships between the flexural rigidity and spacing of the intermediate 
stiffeners and the buckling stress of the stiffened web plate have been obtained. 


One interesting and important feature of these new relationships is that 
they define more clearly than hitherto the difference in the behaviour of single- 
and double-sided stiffeners. 


1. Introduction 


Plate girders are extensively used in one form or another in all types of civil 
and aeronautical construction and, because of the obvious necessity of obtaining a 
correct design procedure, many investigations have been made to examine the 
behaviour of plate girders when under load. As a result, much is now known about 
the various buckling phenomena associated with plate girders. 


However, one very important problem, which has not attracted as much 
attention as others, is the determination of the relationships which exist between 
the size and spacing of intermediate vertical stiffeners and the shear buckling stress 
of the stiffened web plate. The existing experimental evidence, while inconclusive 
in itself, has shown that in practical structures the flexural rigidity of the inter- 
mediate stiffeners has to be much greater than that required by the theoretical 
analyses. This is not surprising when it is appreciated that, to obtain the theoretical 
solutions, it has been necessary to make certain simplifying assumptions which are 
not strictly correct when practical structures are considered. 


This paper presents the results obtained from an experimental investigation 
made in an attempt to obtain the relationships between the size and spacing of 
intermediate vertical stiffeners and the buckling stress of the stiffened web plate 
when it is subjected to shear. Tests were made on more than two hundred different 
plate-stiffener combinations and from the analysis of the results it has been possible 
to derive new empirical relationships. 


Received September 1955. 
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NOTATION 
t thickness of web plate 
t, thickness of attached stiffener leg 
b spacing of intermediate stiffeners or width of unstiffened web plate 
b, clear web plate distance between stiffeners 
‘d, clear depth of web plate 
a, effective aspect ratio 


=b/d, for unstiffened web plates or web plates reinforced by single- 
sided stiffeners 


=b./d. for web plates reinforced by double-sided stiffeners 
D flexural rigidity of unit width of plate= / {12 (1 
E Young’s modulus 
~Poisson’s ratio 
I moment of inertia of stiffener 
EI flexural rigidity of stiffener 


_ stiffener rigidity 
Y= Db _ rigidity of strip of web plate equal in width to stiffener spacing 


I, moment of inertia of the combination of a single-sided stiffener plus 
strip of web plate equal to stiffener spacing, calculated about the 
centroid of the combination 


Y, limiting value of y 
T shear stress 
T., critical shear stress= KD /(d.?t) 
K critical shear stress coefficient 
K, limiting value of K 


2. Background 
2.1. SUMMARY OF EXISTING THEORETICAL KNOWLEDGE 


Theoretical studies“: ?:*» have established that the buckling stress T,, of plane 
unstiffened rectangular panels can be determined from equation (1), in which K 
is a coefficient, the value of which depends upon the superficial dimensions of the 
panel and the type of edge support. 


Kr? 


d.*t 


(1) 


As Fig. 1 shows, the value of the critical shear stress coefficient K increases 
with decreasing values of the aspect ratio b/d... This explains why intermediate 
vertical stiffeners provide such an effective method of increasing the buckling stress 
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Fig. 1. 


of rectangular panels. In view of the importance of obtaining the correct design 
of intermediate stiffeners, it is not surprising that a number of theoretical investi- 
gations have been made to determine the relationships between the size and spacing 
of intermediate stiffeners and the buckling stress of the stiffened web plate. 


Schmieden, Seydel) and Wang have obtained relationships between the 
size and spacing of relatively weak, but closely spaced, transverse stiffeners and 
the buckling stress of the stiffened web plate. The effect of such stiffeners may 
be considered as stiffening the plate in one direction, but not influencing the 
strength in the direction normal to the stiffeners and the problem reduces to a 
study of the buckling of an orthotropic plate. From the exact results obtained by 
Seydel, certain well known design formulae have been developed”: *. 


In 1936, Timoshenko“’ and Way™ obtained relationships between the non- 
dimensional parameter y(=EI/Db) and the shear stress coefficient K of simply- 
supported finite rectangular plates reinforced by one or two transverse stiffeners. 
They showed that if the transverse stiffeners were relatively flexible they would bend 
with the plate and not restrict the wave pattern, but that, by increasing the flexural 
igidity of the stiffeners, a stage would be reached when they would remain straight, 
thus forming a nodal line. From their solution, Timoshenko and Way calculated 
the minimum flexural rigidity which stiffeners should possess if they are to remain 
straight and supply adjacent panels with a simple pin-joint support; the values they 
obtained are plotted in Fig. 2. 
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Fig. 2. 
Timoshenko and Way’s theoretical relationships for the condition that intermediate stiffeners 
supply a simple pin-edge support. 


In 1949, Stein and Fralich“® produced an excellent paper in which they made 
a thorough study of the behaviour of infinitely long, simply-supported plates rein- 
forced by equally spaced transverse stiffeners. Now, theoretical work on plane 
rectangular panels has shown that buckling loads above the true minimum value 
will be obtained if an incorrect, that is ill-chosen or inadequate, expression is used 
to define the deflected form of the plate. In their paper, Stein and Fralich made 
a detailed study of this problem to ascertain which deflection expression would 
provide the correct, that is minimum, buckling load. The K/y relationship they 
obtained, for one of the aspect ratios considered, is plotted in Fig. 3, and it will be 
noted that points of discontinuity occur on the K/y curves. These points of 
discontinuity denote where changes in the governing buckle pattern occur. This 
work of Stein and Fralich showed that the earlier theoretical work of Timoshenko 
and Way overestimated the value of the critical shear stress coefficient K obtained 
with stiffeners of intermediate stiffness, and underestimated the value of K obtained 
with fairly rigid stiffeners. 
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Fig. 3. 
The theoretical K/y relationships derived by Stein and Fralich for infinitely long plates simply- 
supported at the edges and reinforced by stiffeners 0-5d, apart. 


2.2. SUMMARY OF EXISTING EXPERIMENTAL INVESTIGATIONS 


Throughout the past century, a great number of tests have been carried out on 
plate girders in attempts to obtain a better understanding of the behaviour of the 


| web plates when under load. However, on examining the technical literature, it 


will be found that very few of the investigations were primarily concerned with 
the influence of stiffener dimensions upon the stability of the web plate. 


One notable exception was the investigation of Moore“") who, in 1942, reported 
on his extensive tests on two aluminium alloy girders. Unfortunately, Moore failed 
to obtain values of the experimental buckling loads and had to resort to determining 
test loads corresponding to an arbitrary constant value of web or stiffener deflection. 
One great disadvantage of this procedure, and one which Moore himself recognised, 
was that the use, as a measure of the effectiveness of stiffeners, of a constant value 
of web or stiffener deflection at a load corresponding to the theoretical load for 
panels with all edges 50 per cent. fixed, will of necessity involve different degrees 
of buckling for panels of varying size and different aspect ratios. 


Moore, working on these lines, produced the empirical design formulae 


Y= 
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Design details of the BTG/- series of test girders. 


Now, although a designer who might use this formula would not know the buckling 
load of the plate-stiffener combination, he would know that at loads corresponding 
to the theoretical buckling loads for panels with partially (50 per cent.) clamped 
edges, the web plate and stiffener deflections would be small. 


A year later, Scott and Weber?) reported on tests on stiffened shear and 
compression panels. Like Moore, they experienced difficulty in obtaining experi- 
mental values of the buckling load. As a result of their tests they found that an 
effective stiffener required a flexural rigidity greater than that predicted by 
Timoshenko and Way. 


In the United Kingdom, Sparkes“® has conducted an extensive programme 
of research on plate girders and, in the course of his tests, he also found that 
intermediate stiffeners having a flexural rigidity equal to twice Timoshenko’s 
theoretical value failed to divide the web plate into separate panels. 


From the foregoing it will be noted that, although the previous investigations 
failed to provide experimental relationships between the parameter y and the shear 
stress coefficient K, they did establish the fact that intermediate stiffeners designed 
in accordance with Timoshenko’s theoretical relationships would be unsatisfactory, 
in that they would not provide the necessary support to the plate. 
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3, The Experimental Investigation 
3.1. AIM OF INVESTIGATION 


The testing programme was planned to provide information on the following 
points: — 


(i) The relationships between the critical shear stress coefficient K and the 
non-dimensional parameter y for varying values of the ratio 
(stiffener spacing) /(clear depth of web plate). 


(ii) The difference in the behaviour of single- and double-sided stiffeners, 
which also involved an examination of the influence of the width of the 
attached leg. 


3.2. THE DESIGN OF THE PLATE GIRDERS AND SHEAR PANELS 


The plate girders and shear panels were designed to conform as nearly as 
possible with the following requirements: — 


(a) that the web plates be subjected to pure shear, 


(b) that the buckling stresses of the stiffened web plates be well below the 0-1 
per cent. proof stress of the sheet material, 


(c) that failure of the girder should be due to failure of the web plate. 


Requirements (a) and (c) were satisfied by using relatively heavy flange members 
which ensured that the applied bending stresses were small. By constructing the 
girders and shear panels of H 10 WP* high strength aluminium alloy and using 
web plates having high d,/t ratios, requirement (b) was satisfied. 


The design details of the BTG/- series of test girders are shown in Fig. 4, 
while details of the shear panels are clear from Fig. 5, which shows a specimen 
in the testing machine. 


3.3. THE TESTING EQUIPMENT 


The plate girders and shear panels were tested in an Amsler testing machine, 
the general features of which can be seen in Fig. 6, which shows a girder in 
position. The girder was supported on roller bearings and loaded through a 
similar roller bearing, this ensuring that the flange members could expand or 
contract during a test. 


The lateral deflections of the web plate were measured with accurate dial 
gauges graduated in units of either 0-0005 in. or 0-0001 in. To ensure that the 
dial gauges only recorded the lateral deflections of the web plate, they were attached 
to a special frame which was secured to the girder by three adjustable screws 
situated in the plane of the web plate. (For details, see Fig. 6.) This arrangement 
ensured that there would be no relative movement of the dial gauges with respect 
to the central plane of the girder, should the girder twist under load. During some 
of the tests, additional dial gauges which were not attached to the special frames 
were used. Unfortunately, the results obtained with these gauges were often 
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Fig. 5. 
Shear panel in the testing machine. 


unsatisfactory, since the gauges had recorded some overall movement of the girder 
in addition to the lateral deflection. 


During certain tests, contour plots were taken of the buckled shape of the 
web plate. To obtain these contour plots, a series of gauge bars were held parallel 
to the central plane of the girder and the position of the web plate relative to 
these bars was determined before and during the test, using a dial gauge graduated 
in units of 0-0005 in. 


In addition, electrical resistance strain gauges were used to record the strains 
in the web plate and stiffeners. The gauges used were British Thermostat Self- 
Adhesive gauges having an effective gauge length of 4 in. 
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Fig. 6. 
A test girder shown in position in the testing machine. 


3.4. THe Tests 


The testing programme involved tests on 18 plate girders and 12 shear panels. 
Six of the test girders were of riveted construction, and 12, referred to as the 
BTG/- series of test girders, were of bolted construction. No details are given 
of the tests on the riveted test girders, since large lateral deformations were 
formed in the web plates during construction and consequently the buckling 
phenomenon was not well defined. These tests made it clear that little progress 
would be made towards obtaining the K/y relationships unless the web plates 
of future girders were flat, thus ensuring that a well defined buckling phenomenon 
would be obtained. It was for this reason that the BTG/- series of test girders 
were of bolted construction. 


The first girder of the BTG/- series was very carefully constructed and two 
tests were made on it. The first test, done before intermediate stiffeners were 
fitted, was to determine whether, in accordance with Ref. 15, the web plate was 
rigidly clamped at its edges. The web plate was flat and, as a result, extremely 
good plots of load against lateral deflection and load against strain were obtained. 
For the second test on girder BTG/1, the web plate was divided into three equal 
panels by employing two sets of very stiff double-sided stiffeners. During the 
initial stages of this test, electrical resistance strain gauge and dial gauge readings 
were taken at frequent intervals until it was apparent that the panels had buckled. 
Once again a very well defined buckling phenomenon was obtained, as shown 
in Fig. 7. 
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Load/strain plots obtained from test BTG/1/2. 


Girder BTG/2 was fitted with a weaker set of double-sided stiffeners and a 
test similar to the second test on girder BTG/1 was done on it. Once again a well 
defined buckling load was obtained. These tests therefore showed that the degree 
of planeness of the web plate obtained with the method of construction was 
satisfactory. 


The next step was to determine whether good results could be obtained when 
single-sided stiffeners were used, and tests involving single-sided stiffeners were 
conducted on girder BTG/3. At first a set of relatively flexible single-sided 
stiffeners were fitted and the girder slowly loaded, while dial gauge and strain gauge 
readings were taken at frequent intervals. When it was evident that the web plate 
had buckled the girder was unloaded. The resulting zero load readings were then 
compared with the initial zero load readings to determine whether the web plate 
had, in fact, returned to its original position. This was found to be the case. The 
intermediate stiffeners were then removed and, after fitting a slightly stronger set 
of intermediate stiffeners, the loading and unloading procedure was repeated. From 
the experimental readings obtained from both tests, well defined buckling loads 
were obtained. 


The tests had clearly established the fact that, provided that the web plate 
was initially plane, well defined buckling loads could be obtained with both single- 
and doubled-sided stiffeners. 
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Since it was intended that many tests should be done on each of the remaining 


: girders, it was considered that a check on the influence of repeated buckling upon the 


behaviour of a panel would be of value. A repeated loading test was, therefore, 
conducted on an unstiffened shear panel. (For details of the loading frame, see 
Fig. 5.) 


During this test, the panel was repeatedly loaded up to a load corresponding 
to twice the initial buckling load of the panel, the buckling load of the panel being 
determined during certain of these operations. It was found, as was expected, 


| that the buckling load was unaffected by the loading operations. This was of 


interest since, during certain shear tests on curved panels“, it had been found 
that repeated buckling had lowered the buckling stress of the panels. 


In all the subsequent tests on the test girders and shear panels the following 
procedure was adopted. The weakest set of intermediate stiffeners to be used on 
the given girder would be placed in position and the zero load readings taken. 
The girder was then gradually loaded, dial gauge readings, and in certain tests 
strain gauge readings, being taken at frequent load intervals. When it was apparent 
from these readings that the web plate had been loaded up to at least 30 per cent. 
beyond the buckling load, the girder was unloaded. The resulting dial and strain 
gauge readings were then compared with the initial zero load readings. If the dial 
gauge readings returned to within a few divisions of the initial readings (each 
division was 0-0001 in.), it was considered that the web plate had not been 
permanently deformed during the test and that a further test could be done on it. 
The intermediate stiffeners were then removed, another set of stiffeners fitted to the 
web plate, and the procedure repeated. 


4. The Analysis of the Experimental Results 


Altogether 220 different plate-stiffener combinations were tested, of which 125 
involved single-sided stiffeners. The K/y relationships obtained from most of the 
tests followed a well-defined pattern. However, there were a few cases with widely 
spaced stiffeners when, because of excessive scatter of the results, the relationships 
were not so apparent. With regard to the scatter of experimental results, it will 
be appreciated that if a few freak or poor results are obtained, then a much larger 
number of normal results are required before the correct or real relationships 
become apparent. 


An examination of the results obtained from two sets of tests, one involving 
double-sided stiffeners, the other single-sided stiffeners, will illustrate the nature of 
the K /y relationships. 


Twenty-three tests involving single-sided stiffeners were done on three girders, 
BIG/7, 8 and 11, the stiffener spacing being 43 in. and clear depth of web plate 
12 in. The experimental results have been plotted in Fig. 8; the value of the 
critical shear stress coefficient K increases with increasing values of y, until a certain 
limiting value y,, which in this case is 124, is reached, while for values of y greater 
than 124, the value of K remains approximately constant at the mean value of 43. 
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Results obtained when employing single-sided stiffeners at 4% in. centres. 
The clear depth of web plate was 12 in. 


It was found that the expression defining the relationship between K and y 
over the range 0 = y < 124 was 


K=11:2+6-35 ‘ ; ‘ (3) 


In addition 25 tests involving double-sided stiffeners were done on girders 
BTG/7 and 8 (see Fig. 9). As in the previous case, the value of K increases with 
increasing values of y until y reaches a certain limiting value and, for values of 
y greater than the limiting value y,, K remains constant. 


The K/y relationship over the range 0 <y < 216 is 
K=11:2+7:3 ‘ ‘ (4) 


Since the ratio (clear depth of web plate)/(stiffener spacing) was the same in 
both of these tests, the resulting K/y relationships provide some interesting infor- 
mation. Firstly, the K/y relationships are similar in form (see Fig. 10). Secondly, 
the tests have shown that, for a given stiffener spacing, higher critical shear stresses 
can be obtained when using double-sided stiffeners. This is because, with double- 
sided stiffeners, the effective panel width is approximately equal to the clear web 
distance between stiffeners, whereas with single-sided stiffeners it is closer to the 
actual stiffener spacing. 
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The experimental K/y relationships obtained from tests involving double- and single-sided 
stiffeners at 43 in. centres on web plates having a clear depth of 12 in. See Figs. 8 and 9. 
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With reference to this last point, in order to determine the effect which the 
width of the attached stiffener leg has upon the critical shear stress of plates 
reinforced by single-sided stiffeners, some 28 tests were done using stiffeners with 
varying widths of attached legs. It was found that although, for a given value of 7, 
there was a tendency for the value of the critical shear stress to increase with 
increasing stiffener width, this increase was small and could be neglected. Fig. 10 
shows that the K/y relationships consist of two distinctly different sections. Over 
the first section, the value of K increases rapidly with increasing values of y until 
a certain value of y, termed the limiting value y,, is reached. Over the second 
section, that is for values of y greater than y,, the value of K remains constant at K;. 


It has been seen that, over the first section, the K /y relationship can be defined 
by a cubic expression similar to 


K=K,+A (yy)? 
K=K,, for values of y>y1 
where K =critical shear stress coefficient of the stiffened web plate, 


K,=critical shear stress coefficient of the unstiffened plate. 
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Thus, provided that the values of y, and K, were known, it would be a 
relatively simple procedure to calculate the critical shear stress of a plate stiffened 
by identical stiffeners equally spaced along the plate. 


The values of K;, obtained from both the single- and double-sided stiffener 
tests have been plotted against the corresponding values of the effective aspect ratio 
a, in Fig. 11. 


The relationships between K;, and the effective aspect ratio 2, for single- and 
double-sided stiffeners are given respectively by 


—2 
Single-sided stiffeners: K,=8:0+5-7 (4) (5) 


Double-sided stiffeners: K,=7:0+5-6 (2) (6) 


It is interesting to note that these relationships are similar in form to the 
approximate relationships between K and 2, for plates which are either clamped or 
simply-supported on all four edges“”. 


The values of +; obtained from the tests have been plotted in Fig. 12 and the 
relationships between y, and the effective aspect ratio z, for single- and double-sided 
stiffeners are respectively 


Single-sided stiffeners: y,=21-5 (2) (7) 


—2 
Double-sided stiffeners: y,=27-75 -7°5 (8) 


The two K,/z, relationships are very close to each other, the curve for single- 
sided stiffeners lying just above the curve for double-sided stiffeners. It could be 
argued that the reason for this is that the effective panel width for single-sided 
stiffeners is a little less than the stiffener spacing and that for double-sided stiffeners 
is a little more than the clear web distance between stiffeners. If this is the case, 
then the effect would be to increase the gap between the two curves plotted in 
Fig. 12. It is therefore clear that different design rules should be used to determine 
the flexural rigidity which single- and double-sided stiffeners should possess. 


With reference to Fig. 11, since the two curves are so close to each other, it is 
proposed that for design purposes the relationship for double-sided stiffeners 
(equation (6) ) be used for both types of stiffening. 


It is now generally accepted that, if the thickness of the attached stiffener leg 
is less than the thickness of the web plate, the effectiveness of the stiffener is 
reduced. Now, most of the tests in this investigation involved stiffeners having 
attached legs thicker than the web plate, because it was considered that little 
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progress could be made towards estimating the influence of the ¢,/f ratio until 
more was known about the general behaviour of intermediate stiffeners. Most of 
the stiffeners used had a thickness within the range ¢ to 1-6t. 


Therefore, until more is known about the influence of the ¢,/t ratio, the 
empirical relationships given in equations (5) to (8) should only be used when the 
attached stiffener leg is equal to or greater than the thickness of the web plate. 
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4.1. SUMMARY OF THE EMPIRICAL RELATIONSHIP OBTAINED FROM THE 
INVESTIGATION 


The analysis of the experimental results has shown that, for panels loaded in 
shear, the following empirical laws define the relationships between the critical shear 
stress coefficient K, the non-dimensional parameter y, and the effective aspect 
ratio 

K= K, A (y)', 


where K,,=critical shear stress coefficient of the unstiffened plate, 


K =critical shear stress coefficient of the stiffened plate, being equal to 
K, for values of y > yx. 


isaconstant, 
y=ElI/(Db) 
27:75 (a.)~? — 7:5 for double-sided stiffeners 
21-5 — 7:5 for single-sided stiffeners, 
K,=7-0+5:6(2.)~? for both single- and double-sided stiffeners, 
and a,=b,./d. for double-sided stiffeners, 
b/d. for single-sided stiffeners. 


These relationships are only valid for intermediate stiffeners having an attached 
stiffener leg whose thickness is equal to or greater than the thickness of the 
web plate. 


5. Comparison of Theoretical and Experimental Values 


The experimental results support in principle the theoretical work of Stein and 
Fralich. The form of the experimental K/y relationships (Figs. 8, 9 and 10) is 
similar to that obtained theoretically by Stein and Fralich. There is, however, a 
difference in the magnitude of the values. This is understandable when it is 
appreciated that the theoretical values were derived for the case of a plate simply- 
supported at the edges, whereas the empirical values have been obtained from 
tests on practical girders in which the flange members supply a high degree of 
rotational restraint to the web plate. 


An important difference between the theoretical and experimental work is that 
from the experiments it has been possible to distinguish between the behaviour of 
single- and double-sided stiffeners. 


To provide a direct comparison between the empirical and theoretical relation- 
ships, curves have been obtained by applying equations (4), (6), (7) and (8) to the 
case of an infinitely long plate reinforced by stiffeners spaced at d., 0-5d, and 0:2d.; 
they are plotted in Fig. 13, together with Stein and Fralich’s theoretical curves. The 
degree of agreement increases as the stiffener spacing decreases, this no doubt 
being due to the fact that the increase in the buckling coefficient K obtained with 
the practical girders, due to rotational restraint supplied by the flanges, becomes 
proportionally less with a decrease in stiffener spacing. 
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6. The Design of Intermediate Vertical Stiffeners 


In view of the previous lack of information regarding the influence of the 
flexural rigidity of intermediate stiffeners upon the buckling stress of plate-stiffener 
combinations, it is not surprising to find considerable differences in the various 
design laws in use at the present time (Fig. 14). 


It will be noted that the values of y recommended in the Royal Aeronautical 
Society’s Data Sheet“* 02.03.14, and which are very close to the values recom- 
mended in the German Specification D.4114°, are much smaller than the empirical 
values given by equations (7) and (8). The data sheet does recommend that when 
the edges of the web plate are subject to angular restraint, the value of y should be 
increased by approximately 50 per cent. However, even after making this adjust- 
ment, for values of z. greater than 0-5, the data sheet value is low. Such stiffeners 
will no doubt deflect with the plate and fail to form a nodal line. The same 
criticism also applies to the design law which is used in certain American 
Specifications®: 21), 


In addition to its function of increasing the buckling stress of a plate, an 
intermediate stiffener must help to support the flanges when the web plate is loaded 
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beyond its buckling load. Therefore, intermediate stiffeners on shear panels 
designed to operate at high values of T/T.., have to be strong and rigid in order to 
be capable of carrying the vertical component of the diagonal tension field. This 
requirement has meant that most of the stiffeners used in aircraft structures have 
possessed an inertia greater than that required by Data Sheet 02.03.14. However, 
now that more shear panels are being designed to operate at low values of 7/T.;. 
it has become necessary to know the relationships between the size and spacing of 
intermediate stiffeners and the buckling stress of the stiffened panels. 


In the United States the National Advisory Committee for Aeronautics” 
has produced a design procedure for shear panels based upon the extensive tests 
conducted by Kuhn®?:**) and others®*:*5), This design procedure does not contain 
any clause which states the minimum flexural rigidity which stiffeners should 
possess. There are, however, two charts from which, together with a formula, 
the buckling stress of a stiffened panel can be calculated. Instead of employing 
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the flexural rigidity of the stiffener, the parameter ¢,/t is used and two curves are 
supplied, one for single-sided stiffeners and one for double-sided stiffeners, which 
provide values of a restraint coefficient in terms of the parameter ¢,/t. Using 
this design procedure, the values of the buckling coefficient K plotted in Fig. 15 
can be obtained. Also plotted in Fig. 15 is the empirical relationship given by 
equation (6). 


According to the American design procedure, for given values of the parameter 
t,/t and the aspect ratio «., a much lower value of the critical shear stress coefficient 
will be obtained with single-sided stiffeners than with double-sided stiffeners. The 
investigation conducted by the author has indicated, however, that the value of K 
obtained with single-sided stiffeners will be very close to the value obtained with 
double-sided stiffeners spaced so as to provide the same value of 2,. In fact, in 
Fig. 11 the curve for single-sided stiffeners actually lies just above the curve for 
double-sided stiffeners. 

Another point of interest is that, as a result of the tests, it is doubtful whether 


there will be so large a gain in the value of K by using values of ¢,,/t greater than 1, 
as is indicated in Ref. 15. 


A further series of tests is now being done to investigate more fully the 
parameter ¢,/t and, now that the general relationships between the variables K, y 
and 2, are known, it is expected that it will be possible to determine quite accurately 
the influence of the parameter ¢,/t. 


7. Conclusion 


The empirical relationships between the buckling stress of a stiffened panel 


_ and the flexural rigidity and spacing of the intermediate stiffeners, which have been 


presented in this paper, will enable designers to predict with greater accuracy the 
buckling load of uniformly stiffened plates subjected to shear. One important 


| feature of the new relationships is the different formulae which have been obtained 
_ for single- and double-sided stiffeners, thus defining more clearly than hitherto the 


difference in the effect of the two types of reinforcement. 
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Photoelastic Investigation on Plates with Single 
Interference-Fit Pins with Load Applied to 
Plate Only 


H. T. JESSOP, C. SNELL and G. S. HOLISTER 


(Department of Civil and Municipal Engineering, University College, London) 


SUMMARY: The introduction of an interference-fit pin in a circular hole in a 

flat plate is found to result in a reduction in the stress concentration factor at 

the hole boundary under tension applied to the plate. The S.C.F. decreases 

rapidly with increase in the ratio of interference stress to applied tension, and 

decreases also with increase in the ratio of hole diameter to width of plate. 

Very little difference in the S.C.F. is found for different ratios of Young’s 
modulus of pin and plate. 


I. Introduction 


This is the third stage of the photoelastic tests in a programme of research 
initiated by the Royal Aeronautical Society and outlined in Ref. 1. The first and 


| second stages dealt with the stress distributions in tension plates with unfilled holes, 


and holes filled by a push-fit pin, respectively, and the results were published in 


| Refs. 2 and 3. This investigation deals with tension plates with circular holes filled 


by interference-fit pins, with the load applied to the plate only. 


NOTATION 
width of plate 
nominal diameter of empty hole 
diameter of pin 
Young’s modulus of plate [E (bar)]* 
Young’s modulus of pin [E (pin)]* 
Poisson’s ratio of plate 
modulus of rigidity of plate=E/{2(1+v)} 
P.q principal stresses 


| ‘The symbol in square brackets is the one which has been used in previous papers 
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‘ 


o,,09 radial and circumferential normal stresses 
t mean tensile stress in undrilled plate due to applied tension [r,]* 


Dt/(D—d), nominal mean tensile stress across section through the 
centre, due to applied tension [?’,]* 


tm | maximum tensile stress on hole boundary 
Sm = 4(p—q)m, Maximum shear stress on hole boundary. 


The stress concentration factors derived here are given in terms of the 
increase in peak tensile or shear stress produced by the applied tension. 


II 


Thus, S.C.F. for tensile stresses = d1,,/t 
S.C.F. for shear stresses = 6s,,/f. 


| 


2. Conditions of Tests 
2.1. THE MODELS 


As in previous tests the models were Araldite plates, 1 in. wide and } in. thick, 
with four hole sizes: — 4, }, 3 and 4 in. diameter respectively. 


Most of the tests were made with pins of brass and Araldite, but a few were 
made with a type of filled Bakelite. The ratio of the Young’s moduli for the 
pins were: — 

Brass Bakelite Araldite 
E/E = 30 2°15 1. 


*The symbol] in square brackets is the one which has been used in previous papers. 
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PHOTOELASTICITY 


1.2. DEGREES OF INTERFERENCE 


Two degrees of interference were examined for each size of hole. These were 
approximately 


0:3 per cent and 0-7 per cent for the brass pins 
0-4 per cent and 0-9 per cent for the Araldite and Bakelite pins. 


These values were found to give comparable amounts of interference stress 
with the different pin materials and they are specified as “low” and “high” 
interference respectively. 


These two degrees of interference would produce in light alloy plates maximum 
elastic shear stresses at the hole boundary of approximately 21,000 Ib. /in.* and 
50,000 Ib. /in.? respectively. 


2.3. APPLIED TENSIONS 


In the first tests, two different tensions were applied to each plate for each 
degree of interference. In relation to the strains produced these would correspond 
to mean elastic tensile stresses in light alloy plates of 15,000 lb./in.? and 
30,000 Ib./in.? in the undrilled plate. 


Later, in order to trace the changes in S.C.F’s more thoroughly, it was found 
necessary to subject the plates to a range corresponding to 0 to 40,000 Ib./in.” in 


| the light alloy. 


| 24. METHOD OF TESTS 


The fringe patterns under the various loads were photographed in light and 
dark fields on 35 mm. film, and the boundary shear stresses were obtained by 


extrapolation from graphs plotted from enlargements of the photographs. 


To obtain the tensile stresses on the boundary, sets of isoclinic fringes were 


traced in the projecting polariscope and the separate stresses were evaluated by 
| graphical integration of the equations of equilibrium along the horizontal and 
| vertical axes of symmetry. 


3. Stresses due to Interference Fit of Pin 
3.1. DEGREE OF INTERFERENCE 


Preliminary tests showed that measurement of the diameters of pin and hole 
failed to give a reliable estimate of the interference stresses to be expected in a 
plate, and that this was due to inability to measure the hole diameter to a sufficient 
degree of accuracy. Further tests with a series of pins measured to an accuracy of 
0:0001 in. and inserted in the same hole showed that the shear stress on the hole 
boundary varied linearly with the diameter of the pin. Extrapolation from a graph 
then yielded the point at which the production of zero stress indicated the effective 
hole diameter. In every case tested this was found to be somewhat greater than 
the measured hole diameter. 
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In the tests described here therefore, the holes and pins were turned to give as 
nearly as possible the required degree of interference stress, and this stress, as 
measured from the photoelastic fringe pattern, was taken to indicate the degree of 
interference obtained. 


3.2. THE SHEAR STRESSES 


Figures 2 and 3 show the photoelastic fringe patterns in a plate with a } in. 
diameter hole, with low and high interference respectively. The dark fringes are 
loci of equal shear stress, the stress at each fringe being proportional to the fringe 
number as shown by the numbers on the photographs. 


In all cases the stress around the hole boundary was subject to small variations, 
through lack of perfect circularity of pin and hole. Also the presence of the pin 
tended to mask observation of the boundary itself, so that the value of the true 
boundary stress had to be found by extrapolation from a graph of the stresses along 
a radial line. The values of the peak fringe numbers so obtained are liable to 
errors of the order of one-fifth of a fringe (i.e. about 5 per cent in the most 
unfavourable case). Within this limit of error it was found that :— 


(i) The value of the maximum shear stress was constant around the hole 
boundary. This applied for all hole sizes and for both low and 
high interference. 


(ii) The shear stress on the hole boundary for a given degree of interference 
was the same for all hole sizes. 


(iii) The magnitudes of the maximum shear stresses agreed with those given 
by the theoretical solutions for an oversized pin in an infinite plate. 
(See Appendix). 


3.3. THE SEPARATE STRESSES 


(i) At the hole boundary on the vertical axis of symmetry (points C, C’ in 
Fig. 1) the circumferential tension and radial compression are equal in 
magnitude, i.e. 7,=\—o, as in the infinite plate solution. 


(ii) At the hole boundary on the horizontal axis, (points A, A’ in Fig. 1), the 
circumferential tension is greater and the radial compression less, the 
difference being small for the 4 in. diameter hole, and increasing with 
the hole size. 


The increase in «, and the decrease in | 7, | due to finite width of the 
plate were, in fact, equal and so caused no variation in shear stress 
around the hole. 


4. Qualitative Effects of Application of Tension to the Plate 

Figures 4 and 5 show photoelastic fringe patterns for the plate with a } in. 
diameter hole under tension, Fig. 4 when a small tension is applied with the high 
interference pin, Fig. 5 when a large tension is applied with the low 
interference pin. 
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Fig. 2. Fig. 3. 
Fringe pattern due to interference fit of Fringe pattern due to interference fit of 
pin. 3 in. diameter hole, brass pin, low pin. ; in. diameter hole, brass pin, high 


interference. interference. 


rterly 
Fig. 4. Fig. 5. 
Fringe pattern in plate with high Fringe pattern in plate with low 
interference and low applied tension. interference and high applied tension. 
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(pane TENSILE STRESS 
/ IN UNDRILLED PLATE ) 


SHEAR STRESS (UNITS) 


Fig. 6. 


Shear stresses on horizontal section AB in piate with } in. hole with low interference under 
zero, low and high loads. 


The two points of zero shear stress above and below the hole in Fig. 4 are 
approximately the points at which the applied vertical tension has neutralised the 
vertical compression due to interference. In the region between these points and 
the hole the plate still retains, in this case, some of its initial circumferential tension 
and radial compression. 


As the ratio of applied tension to interference stress increases, these points 
move closer to the hole, the pressure between pin and plate at the top and bottom 
of the hole decreases, and slip begins to occur. Eventually the plate begins to 
“lift” from the pin, and in these top and bottom regions the hole boundary 
becomes “ free” and is subject only to a small circumferential compressive stress. 
This is the condition illustrated in Fig. 5. 


The foregoing sequence is followed as the applied tension increases, for both 
low and high degrees of interference, but in this series of tests the tensions applied 
did not reach a sufficiently high value to produce lift with the high interference pin. 


5. The Shear Stresses due to Applied Tension 


5.1. THe Errect or “ Lirt” 

With the low interference pins, lift first became visible under an applied 
tension ¢ approximately equal to the initial circumferential tension produced by the 
interference. Lift at top and bottom points must have begun at a somewhat 
lower load. 
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HIGH 
INTERFERENCE 


Low 
INTERFERENCE 


SHEAR STRESS (UNITS) 


MEAN APPLIED TENSION ACROSS SECTION THROUGH HOLE 
2 3 4 

Fig. 7. 
Shear stresses on the hole boundary at A under low and high loads. 


In the early tests, in which two loads only were applied to each plate, a 
marked increase in S.C.F’s on the horizontal section was observed under loads 
which produced lift. The later tests under progressively increasing load showed, 
however, that the changes in S.C.F. were continuous and that there was no sudden 
change in either the S.C.F. or its rate of increase corresponding to the actual 
occurrence of lift. 


5.2. THE SHEAR STRESSES ON THE HORIZONTAL SECTION 


(i) Fig. 6 shows the distribution of shear stress* on the horizontal axis of 
symmetry for the low degree of interference under zero, low and high loads 
for the 3 in. diameter hole. The distributions are qualitatively very 
similar for the other hole sizes. 


Table I and Fig. 7 show the shear stresses* at the point A under both 
loads for all hole sizes. Also shown, in broken lines on Fig. 7, are the 
values for the } in. and 3 in. holes under gradually increasing load. 
(These latter curves were obtained from different models from those used 
in the low and high load observations.) 


*The unit of stress corresponds to 15,000 Ib./in.? in a light alloy bar. 
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Fig. 8. 
Shear stress concentration factors at A for low and high loads. 


TABLE I 


SHEAR STRESS AT HOLE BOUNDARY ON HORIZONTAL SECTION 


Shear stress at A 


(units)* (units) 


Hole | 
diameter Low High i Low 
(in.) load load load 


| Increments 
| 
| 


(i) Low interference 


| 
| 


215 
215 
2°15 
2°10 


3°85 
3-40 
3-25 
3°85 


3-65 | 0°70 


3°65 
3°45 
3°85 


0°85 
0-70 
0-60 


(ii) High interference 


4-40 
4-05 
2°85 
4°50 


| 
| 


0°50 
0-15 
0°25 
0°30 


0-60 
0°65 
0°45 
0°30 


0°45 
0°10 
0-15 
0-15 


0°45 
0°30 
0:25 
0-25 


*In this and subsequent tables the unit of stress corresponds to 15,000 Ib./in.? in a light 


alloy bar. 
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TABLE II 


SHEAR STRESS AT HOLE BOUNDARY ON VERTICAL SECTION 


Shear stress at C Increments 
(units)* (units) 


Hole 
diameter Low High Low High 
(in.) © load load load load 


(i) Low interference 
0°66 — 0:05 —0°75 —1°50 
0°60 0:06 —0°70 — 1:25 
0°66 0:26 —0°80 —1:20 
0°60 0°33 —0:90 —1:15 


(ii) High interference 
2°56 1:85 —0°75 — 1-50 
2°30 1-46 —1:10 —1-80 
235 1:46 —0:90 —1:75 
2°45 1:60 — 1:00 —1:85 


On Fig. 7 are also shown the shear stress graphs for a plate with push-fit 
pin. For the latter stages of loading, with the low interference pin, i.e. 
after lift has occurred, the rate of increase of stress with load approaches 
that in the push-fit case. 


Fig. 8 shows the values of S.C.F. as given in Table I, together with 
those for the push-fit pin and for the theoretical solution for a rigid pin 
“fixed ” in an infinite plate. (See Appendix). 


5.3. THE SHEAR STRESSES ON THE VERTICAL SECTION 


Table II shows the shear stresses on the vertical section at the point C on the 
hole boundary for the same series of tests as before. The addition of tension 
causes a reduction in the shear stress throughout the range of loading applied. 


While there is a general tendency for the S.C.F. to decrease slightly numerically 
with increase in hole diameter, this effect is small and irregular, and is probably 
influenced by small differences in the contact conditions with different pins. 


There is no significant difference in S.C.F. for either change of load or change 
of degree of interference until after the occurrence of lift, when a small reduction in 
the factors occurs with increasing load. 


*The unit of stress corresponds to 15,000 lb./in.2 in a light alloy bar. 
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TABLE III 


TENSION AT HOLE BOUNDARY ON HORIZONTAL SECTION 


(units)* (units) 


Hole 
diameter Zero Low High Low High 
(in.) load load load load load 


Tensile stress at A | Increments 


(i) Low interference 
1°57 2:13 4:27 0°56 2:70 
1-43 2°54 4°54 3°11 
1-70 2°47 4-07 0°77 
2°20 2°60 4-80 0°40 


(ii) High interference 
3°56 5-20 0°84 
3°60 4-60 0:20 
3°60 3-87 0:07 
4-80 5-60 | 0°14 


The Tensile Stresses due to Applied Tension 


6.1. THE TENSILE STRESSES ON THE HORIZONTAL SECTION 


These are shown in Table III. The values obtained for stress increases under 
load are subject in this case to the errors involved in two graphical integrations— 
one for the tension, due to interference, and one for the total tension under load 
—in addition to actual observation errors. The values of S.C.F., therefore, show 
more irregularity than those for the shear stresses; but they show the same general 
trends as the shear stresses: —a general decrease with hole diameter, as would be 
expected, and a marked increase when the ratio of applied load to initial interference 
stress becomes high. 


In the most unfavourable case, that of the 4 in. hole with approximately 0-3 per 
cent interference and high load, the S.C.F. is still only about one half of that for 
the same hole with a push-fit pin. 


6.2. CIRCUMFERENTIAL STRESSES 


Table IV gives the results for the circumferential stress at points C, C’ on the 
vertical section. As for the shear stress at the same points, the S.C.F. is consistently 
negative. Under high loads with low degree of interference the circumferential 
Stress may change sign and become a small compression. Accurate evaluation of 
the separate stresses becomes difficult in this case, however, owing to the difficulty 
of integrating through the isotropic points, and the values given can only be taken 
as indicating the order of the tensile S.C.F. at these points. 


*The unit of stress corresponds to 15,000 Ib./in.? in a light alloy bar. 
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TABLE IV 
TENSION AT HOLE BOUNDARY ON VERTICAL SECTION 


Tensile stress at C Increments 
(units)* (units) 


Hole 
diameter Low High Low High 
(in.) | load load load load 


(i) Low interference 
1:20 —0:10 —0°37 —1°50 
0:93 0°13 —0°37 —1°70 
17 0°53 —0:30 —0:90 
1:90 0°67 —0-60 —0°83 


(ii) High interference 
2°94 2:00 —0°40 — 1-34 
2°53: —0°74 —1:27 
3-06 —0°18 —0°64 
2:93 —0°54 —1-40 


7. Comparison of Results for Pins of Different Materials 


Plates with pins of brass and Araldite were tested under progressively 
increasing load with both degrees of interference for three hole sizes, 4 in., 2 in. 
and 4 in. diameter. 


The plate with 4 in. diameter hole was also tested with a Bakelite pin. 


Figures 9 and 10 show the shear stresses at A and C respectively for the 3 in. 
hole. The other hole sizes gave very similar graphs. 


Figure 11 shows graphs of the S.C.F’s at the point A for all the tests. These 
values were derived from graphs similar to those of Fig. 9, and represent in each 
case the ratio of total change in shear stress to applied tension. 


In all except two cases the difference between the S.C.F’s for brass and 
Araldite pins is nowhere greater than 0-2. (The two exceptions occur on the 
vertical section in the early stages of loading with high interference, and may well 
be attributable to some sort of “settling down” process at the contact surface.) 


8. Summary of Results for All Tests 


The two critical points on the hole boundary are those at the extremities of the 
vertical and horizontal diameters. In the former the S.C.F. is always negative, so 
that the maximum stresses there, both shear and tensile, are those due to the inter- 
ference fit of the pin. Also the S.C.F’s at these points are only substantially greater 
than those on the horizontal diameter in the case of small loads. On the horizontal 
diameter the interference stresses are the minimum stresses, and under medium or 


*The unit of stress corresponds to 15,000 1b./in.? in a light alloy bar. 
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Fig. 9. 
Shear stresses at A under progressive loading for brass and Araldite pins. 
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Fig. 10. 
Shear stresses at C under progressive loading for brass and Araldite pins. 
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Fig. 11. 


Stress concentration factors at A for all load and interference conditions. 


308 The Aeronautical Quarterly 


SHE 
Sc 
8s 
Q 2 4 5 
h 
d 
d 
if 
a 
7 | 


rterly 


PHOTOELASTICITY 


IN HOLE 


3 
6 IN.HOLE HIGH 
INTERFERENCELINT. 


"4 SIN.HOLE 
° 
‘4 inf 


BRASS, ARALDITE AND 
BAKELITE PINS 


© 


2 [Sm 


Fig. 12. 


Shear stress concentration factor at A plotted against ratio of interference stress 
to applied tension. 


high loads the S.C.F’s are of the same order of magnitude as those on the vertical 
diameter. Attention was therefore concentrated upon the points on the horizontal 
diameter where, in addition, experience had shown that fatigue failure almost 
invariably occurred. 


The trend of the results obtained suggested that the S.C.F. was related, at 
least qualitatively, to the ratio of degree of interference to magnitude of 
applied load. 


Values of S.C.F. on the horizontal section were, therefore, plotted against the 
ratio of interference stress to mean applied stress across the section through the 
hole, and these graphs are shown in Figs. 12 and 13. In Fig. 12 the points plotted 
include those for all three pin materials, but values of the tensile S.C.F. were 
obtained only in the brass pin cases. 


The accuracy of the S.C.F. as estimated from possible experimental errors is 
of the order of + 0-1, so that the curves drawn in these graphs give a very fair 
representation of the relation as deduced from the photoelastic tests. 


The facts that the values of low and high degrees of interference conform so 
well to the same curves, and that these curves tend toward the points for zero 
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Tensile stress concentration factor at A plotted against ratio of interference stress 
to applied tension. 


interference as obtained from the push-fit pin tests’, indicates that the ratio of 
interference stress to applied tension is in fact a suitable criterion to adopt for 
estimation of S.C.F’s in this problem. 


The significance of these results in the case where the applied tension alternates 
between zero and a known value r is probably best seen by plotting the “mean stress- 
level” (i.e. the mean of initial and final values of stress), and the amplitude of the 
oscillation about this mean value. These will be given by (x+4y)f and Gy)? 
respectively from Figs. 12 and 13. In Fig. 14 this has been done for the shear stress 
for the 3 in. diameter hole. 


It is obvious that the stress level rises very rapidly at high degrees of 
interference and that the useful part of the graph lies in the region of comparatively 
low ratio of interference to applied load. 


9. Effect of Cyclic Application of Loads 


Figure 15 shows the shear stresses on the horizontal hole boundary under 
sequences of loading and unloading. The general effect is to produce a flattening 
of the stress-load curve in the region of low applied load. This is more pronounced 
with the high degree of interference. The stresses under the higher loads appear 
to be unaffected in both cases. 
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Mean value and amplitude of oscillations of shear stress at A for load varying from zero 
to ¢’, plotted against degree of interference. 


The effect of repeated loading and unloading would be, therefore, to raise the 
mean stress level slightly, while reducing the amplitude of the oscillation. 


It should be noted, however, that these effects are likely to vary with different 
surface contact conditions between pin and plate, and the results cannot be taken 
as applying quantitatively to pins and plates of different materials. 


10. Effect of Plastic Yield in Prototype 


The results in these photoelastic tests are valid only so long as the stress 
distribution in the plate is an elastic one. ; 


If in the prototype the interference is so high that an annulus of the plate 
surrounding the hole has undergone complete plastic yield, two modifications will 
be introduced : — 


(i) The amount of vertical tension which the material above and below the 
hole can sustain without approaching the conditions in which “ slip” and 
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Shear stresses at A under cycles of loading and unloading, with low interference. 


“lift” occur is limited to the maximum radial compressive stress in the 
plastic annulus. For loads which produce tensions above this limit it is 
to be expected that a greater portion of the total tension will act across 
the horizontal section of the plate than in the corresponding elastic 
distribution. 


The material in the “ plastic” annulus at the sides of the hole cannot 
sustain as much additional circumferential tension, so that, whatever the 
magnitude of the applied load, more of the additional stress must be borne 
by the remainder of the horizontal section and the area of plastic yield 
will be extended. 


Qualitatively it would appear that plastic yield in the plate, whether 
produced by the interference alone or by the total effect of interference 
and load, would tend to reduce the fatigue strength. 
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Appendix 


Some of the results of the tests have been compared with the following theoretical 
solutions due to A. C. Stevenson. 


(a) Theoretical Solutions for Pins with Interference Fit in a Circular Hole in an 
Infinite Plate. 


(i) For a rigid pin of diameter (d+ ¢) inserted in a hole of diameter d in an 
infinite plate, the stresses at the hole boundary are given by 


where G =the rigidity of the plate= — = 
ere G= igidity of the p °=7i4y° 


(ii) For a pin of the same material as the plate 


(iii) For a pin of modulus E, and Poisson’s ratio v, in a plate of moduli E, vy, 


For a steel pin (E, =30 x 10°, v, =0-25), in a light alloy plate (E=10 x 10°, 
yv=0-25), this becomes 
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(b) Theoretical Solution for a “ Rigid” Circular Inclusion “ Fixed” in an Infinite 
Plate under Tension T. 


(This solution assumes no displacements of the plate at the hole boundary). 


At the edge of the hole on the horizontal section, 


At the edge of the hole on the vertical section, 


> 2 


For an Araldite plate (v—0-38), these become : — 


At the horizontal section 
og= —0-°13T, o,= —0-22T, 
At the vertical section 
o .=1-46T, = —0-285T. 


The stresses at the horizontal section in a plate with a very high interference-fit 
pin might be expected to approach the values given by this solution, but the S.C.F.’s 
as given in Tables I-IV show that, even for the case of high interference and low load, 
the stresses, although small, are still considerably above these values. 


On the vertical section the theoretical solution gives a positive S.C.F. for the tensile 
stress ow», indicating that the boundary conditions there with the interference-fit pin are 
very different from the state of “ fixation” assumed in the theory. 
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An Energy Approach to Climb Performance 


Estimation of a Turbo-Jet Transport Aircraft 


RAMCHANDRA B. TAMBOLI 
(de Havilland Aircraft, Hatfield*) 


SumMMaRY : The climb performance of a turbo-jet transport aircraft is considered, 
using the concept of the energy height. The analytical formulae used in 
evaluating the time needed, distance covered and fuel consumed during climb 
under various conditions are developed. It is concluded that the application of 
this energy approach method to the climb performance estimation of a turbo- 
jet transport aircraft will prove satisfactory from the points of view both of the 
accuracy of the results and of saving some of the considerable time needed for 
the climb calculation by the conventional method. A numerical example is 
given to illustrate the usefulness, applicability and flexibility of this method. 


1. Introduction 


The operating data to be supplied with a transport aircraft in connection with 
the climb performance includes information on the fuel used, distance covered and 
time taken during the climbs to operating altitude. This information will be 
required to cover the appropriate ranges of those parameters which cause 
appreciable change in the quantities. 


With a turbo-jet transport aircraft atmospheric temperature is of importance 
as well as take-off weight and it may be necessary to tabulate information in 
perhaps twenty to thirty different combinations of weight and temperature. The 
time taken to evaluate the climb performance under these various conditions by 
a conventional method (for example, Royal Aeronautical Society Performance Data 
Sheet EG3/1) is considerable. 


The method described here was developed to expedite the climb calculation. It 
has also shown promise of application in other aspects of performance estimation, 
in particular the rapid evaluation of the effects on climb of the changes in any of 
the controlling parameters. It was found convenient to use the concept of the 
energy height throughout the work, thus following the practice established in 
Refs. 1, 2 and 3. 


*Now with de Havilland Aircraft of Canada. 


Originally received February 1955. 
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a, a; bY, by; constants as 
66.4 defined in the 
A, A’, BBC; GR) 

total span of wing (ft.) 


induced drag factor in Cp=Cp,+ C,?/(tAwe) 
(also referred to as aeroplane efficiency factor) 


acceleration due to gravity (ft./sec.?) 


average slope of the curve of climb speed (T.A.S.) against energy height 
(ft. / sec. /ft.) 


average slope of the curve of climb speed (T.A.S.) against temperature 
(ft. /sec. / °C.) 


average slope of the curve of drag against energy height (Ib. /ft.) 
average slope of the curve of drag against temperature (Ib. / °C.) 
average slope of the curve of net thrust against energy height (Ib. /ft.) 
average slope of the curve of net thrust against temperature (Ib. / °C.) 


average slope of the curve of fuel flow against energy height (Ib. / sec. /ft.) 


average slope of the curve of fuel flow against temperature (Ib./sec./ °C.) 


average slope of the curve of weight of aircraft during climb against 
energy height (Ib. /ft.) 


average slope of the curve of weight of aircraft during climb against 
temperature (Ib. / °C.) 


total time to climb to an energy height (sec.) 
aspect ratio (b,?/Sw) 

parasite drag coefficient of aircraft 

lift coefficient 

drag of aircraft (Ib.) 

fuel flow (lb./sec.) 

geometric height of aircraft (ft.) 

energy height of aircraft=H+V*/(2g) (ft.) 
datum energy height from which climb is commenced (ft.) 
relative energy height with respect to H., (ft.) 
pressure altitude (ft.) 

total fuel consumed in climb (Ib.) 

total distance covered in climb (ft.) 

wing area (gross) (ft.”) 


316 The Aeronautical Quarterly 


by 
m 
m 
n 2 
t 
Ay 
Cp: 


CLIMB PERFORMANCE 


temperature in excess of or below I.S.A. (or an arbitrary reference 
temperature selected) (°C.) 


climb speed (ft./sec.) 
weight of the aircraft (Ib.) 
net thrust delivered by the power plant of the aircraft (Ib.) 
angle of climb with reference to ground (radians) 
air density (slugs/ft.*) 
Suffix “0” refers to conditions at the commencement of the climb 


Abbreviations 


E.A.S. equivalent airspeed 
T.A.S. true airspeed 
international standard atmosphere (I-C.A.N.) 


2. The Fundamental Climb Equation 
The fundamental climb equation may be obtained by writing down the 


equation of longitudinal motion of an aircraft climbing at an angle of y in still air as 


(1) 


which may be written 


2 


The energy height is defined as H,=H+ i 


which, after differentiation with respect to time, becomes 


dH, dH V dV 
dt. dt g dt” 


Substituting equation (4) in equation (2) gives 


dH, Vv 


Equation (5) is referred to in this paper as the energy form of the basic climb 
equation. 

The application of this approach to the estimation of climb performance of a 
turbo-jet transport aircraft is described here in two parts, according to whether the 
variation of climb speed, thrust, drag, fuel flow and weight during climb can or 
cannot be assumed to be a linear function of energy height. 
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PART I 
3. Leading Assumption 
In general, for a particular aircraft and at a specified engine r.p.m., the climb 
speed (T.A.S.), thrust, drag, fuel flow and weight during climb can be expressed as 
some function of energy height. In Part I their variations with energy height are 
assumed to be linear or very close to linear. 


In addition, the following particulars of the aircraft are assumed to be known:— 
(i) Drag polar (drag characteristics). 

(ii) Engine r.p.m. for climb condition. 

(iii) Climb speed. (This may be obtained from the analysis of the optimum 
climb technique. In certain cases it may be dictated by the consideration 
of passenger comfort.) 

Engine r.p.m. and climb speed are assumed to be independent of climb weight. 


4. Method 


The thrust and fuel flow under different temperature conditions for which the 
climb calculations are required can be evaluated and plotted against energy height 
and their variations with energy height determined. 


The drag of the aircraft is a function of climb speed and climb weight, of which 
the former is known. The weight of the aircraft during climb is a function of the 
fuel flow and the rate of climb. The rate of climb, in its turn, depends on the weight 


of the aircraft during climb: this is a vicious circle. It is therefore suggested that. 
by using the formula given in Appendix I, first of all a variation of weight of the 
aircraft during climb with energy height should be established by calculating an 
average slope of the climb weight-energy height curve. It is then possible to 
establish the relationship of the drag of the aircraft with energy height based on 
varying weight in climb. 


In short, for a particular aircraft it is possible to establish the relationships 
of the climb speed (T.A.S.), thrust, fuel flow, drag and weight during climb with 
energy height for different climb weights and temperature conditions. 


5. General Case 
5.1. THE Basic EQUATIONS 


Assuming that the thrust, fuel flow, climb speed, drag and weight during 
climb vary linearly with energy height and temperature : — 


Thrust, AT . 6) 
Fuel flow. F=F,+qHx+qAT ‘ (7) 
Climb speed (T.A.S.), V=V,+hHent+hW AT 

Drag, D=D,+mH.n+ mAT 

Climb weight , W=W,+rH.»+rOT ‘ 10) 
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In these expressions X,, F,, V,, D, and W, are the quantities at an arbitrary datum 
energy height, H.., from which the climb is measured; n, qg, h, m, r and n’, q’, h’, m’,r’ 
are the slopes of the thrust, fuel flow, climb speed, drag and climb weight curves with 
energy height and with temperature respectively. AT is the temperature in °C. in 
excess of or below that of the International Standard Atmosphere (I.S.A.), (or an 
arbitrary reference temperature considered). For convenience the concept of the 
relative energy height, Hx, is used throughout this paper. Thus H-» is the energy 
height with reference to an arbitrary datum energy height, H.,, from which the 
climb is commenced (H.n=H.— 


Using equations (6) to (10), equation (5) becomes 


dH ex, 


_ + AT) —(Dy + + m’AT)) (V+ + WAT) 
W,+rHen+rAT 
+ 


W,+rHen+ VAT 
where a’ =h(n—m) 
b’=h(X,—D,)+(n—m) (V,+ WAT) +hAT (n’ 
c’=[(X,—D,)+ AT (n’— (V+ 


5.2. TIME INTEGRAL 


The total time ¢ required to climb from an arbitrary datum energy height, H.,, 
to any other energy height H. is given by (using the concept of the relative energy 
height already explained) 


Hon 
+ 


0 


The solution of this integral depends on whether 
-(i) b? >4dc’ 
(ii) b? <4ad’c’ 
or (iii) b? =4dc’. 
A study of equation (11) will lead to the following conclusions : — 
constant c’ is always positive (for climb X, > D,) 
constant a’ is negative as n is negative and usually greater than m in magnitude 
while h is positive 
constant b’ may be positive or negative. 


Cherefore b’* > 4a’c’ (in practical cases). 
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Hence the time integral becomes 


W,+rAT — rb’ /(2d’) ) 
/ log. (14 bg + log. 


(14 
r b’ 
+ 2a’ log. (1 Her 


where d’= /(b” —4a’c’). 


5.3. DISTANCE INTEGRAL 


The total distance S covered during climb from an arbitrary datum energy 
height H., to any other energy height H. is given by 


Hop 


aad + PAT 


0 


_ (W,+rOT k’ 
A Jog. (1+ I’ Hex rf Aer 


where k’=n-—m 
l'=(X, —D,)+ AT (n’—m’). 


5.4. FUEL INTEGRAL 


Likewise the total fuel Q consumed in climb from an arbitrary energy height 
H., to any other energy height H, is given by 


Q- +rHen+r QT) 


, 9 , , dH. 4 


0 
Hop 
or Q= [4 2a’ 2a” + 


= E’b’ E’ 


+- 


A’=F W,+ (QW, +F OT) 
(q’r+rq) 
E’=rq. 
It may be noted that the integral part and the logarithmic term of equation (16) are 


already evaluated in the time integral; advantage should be taken of this while 
performing a numerical calculation for the fuel integral. 
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6. Special Cases 


6.1. Case I (A): AT STANDARD TEMPERATURE (AT =0) 
Under this condition equations (6) to (10) become 


X=X,+nHer 

F=F,+qHer 

V=V,+hA ex 

D=D,+ mH... 

W=W,+rH.-x. 
Energy climb equation (5) becomes 


D)V/W 


{(X, — (Dy + mH. »)} (V,+hH-r) 
W,+rHer 


_ + bH +c 
W, + 


where a=a’=h(n-™m) 
b=h(X, —D,)+V,(n—m) 
c=V, (X,—D,). 


The time needed, distance covered and fuel used in climb from an arbitrary 
energy height H., to any other energy height H, are evaluated as in the general 
case already discussed. The final results are 


Time Integral 


Hx ) + log, 


1 
b 


d= /(b? —4ac). 
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Distance Integral 


rl log.(1 : Hon) + Bn (24) 
where k=k’=n-—m 
l= X,-D,. 
Fuel Integral 
Hor 


0 


W.q+F,r ( b a rq 
( log. 1+ Hon + Hen a Hex (25) 
Once again the integral part and the logarithmic term of equation (25) will have 
been evaluated in the time integral; advantage should be taken of this while 
performing the numerical calculation for the fuel integral. 


6.2. Case I (B): ConsTANnT CLIMB WEIGHT (dW /dH.,.=r=0) 


Further, if the weight of an aircraft during climb is assumed constant (and 
equal to either the weight at the commencement of climb or an average weight 
during climb), that is, dW /dH.,=r=0, the time, distance and fuel integrals become 
as follows :— 


Time Integral, from equation (23), 


2 1 
W, log. (4 = Hx + log. (26) 
§ b+d 
Distance Integral, from equation (24), 
W,, k 
$= Jog. (1+ Hex). ; 


Fuel integral, from equation (25), 


Hep 
dH ex 


_W _ dil W.q ( b a 
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6.3. Case II 
A particular case with the following assumptions is now considered : — 
(i) standard reference temperature (AT =0), 
(ii) climb speed : constant E.A.S. 
(iii) constant engine r.p.m. 
(iv) climb weight constant during climb (see note below). 


By virtue of these assumptions the forms of equations (17) to (19) remain 
unaffected, while equations (20) and (21) become 


D=D, (constant, neglecting Reynolds number and compressibility effects) 
W = W, (constant). 


Energy climb equation (5) becomes 


dHex + — Vo + 


=(X-D)V/W= W, 
— +b Hen+ 
W, 
where a,=nh 
b,=nV,+h(X,—D,) 
c,=c=(X,—D,) V,. 
Time Integral 
Her 
W, 
a Hex’ +b Aer tc, dH ex 
Wo! log (i+ ) +log 29 
( 
where d,= V(b,’ —4a,c,). 
Distance Integral 
Hop 
W, 
S= 
§= log. 1+ X,-D, ). ‘ (30) 
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Fuel Integral If 


H 


W, (F,+qHex) 


L 
H, 
_ Wogb, F ( b, a, 
(FW, 2a, a,Hen? +b 2a, log. {1+ ent C, Hex 
(31) 
The integral part of equation (31) will have been evaluated in the time integral. 
Note: If a variation of climb weight with energy height is assumed, then the F 
drag, which is a function of weight, will also vary with energy height. Case II 
then becomes identical with Case I (A) already discussed. 
6.4. Case III F 
Another particular case with the following assumptions is considered : — 
(i) standard reference temperature (AT =0), 
(ii) constant engine r.p.m. I 
(iii) climb speed : constant T.A.S. 
(iv) climb weight is assumed to vary with energy height. (Expressions for 
constant weight will be derived as corollaries). 
By virtue of these assumptions the forms of equations (17), (18), (20) and (21) I 
remain unaffected while equation (19) becomes 
V=V,, (constant). 
Energy climb equation (5) becomes ( 
+ Hex) (Do + V 
=(X-D)V/W= W.+1Hex 
Time Integral 


Hor 


_{_ +rHe 
W,(n—m)-r(X,-—D,) 
(n— my? log. [ 1+ Ha] + Viin-m == (32) 
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If the weight is assumed constant during climb (r=0), then 


W, [ _n-m 
t= log. | 1+ Hp |. : (32a) 


0 0 


Distance Integral 


) Hep 


0 
S=V, [time integral (32) ]. ‘ 
a For constant weight during climb (r=0), 
S=V, [time integral (32a) ]. ‘ P . (33a) 
Fuel Integral 
Hep 
V, (X,—D,)+(n—m) Hea} 
For constant weight during climb, 
H, 
for = | oe (F, +qHex) W, 
Q {(X, —D,)+ (n—- m) Ax} (Aa) 
0 
(21) Both equations (34) and (34a) can be evaluated in the usual way. 
PART II 
In those cases where the main assumptions in Part I, namely, that the variations 
of thrust, climb speed, fuel flow and drag are linear with energy height, are found 
hard to justify, the following approach is suggested. 
7. Method 
Where the climb calculations are required for various climb weights and 
temperature conditions by using the method described in Part II, it is advisable to 
select a reference climb weight and temperature such that they lie roughly in the 
middle of the range of the climb weights and temperatures for which the calculations 
are required. For example, if it is expected to perform the climb calculations for 
a climb weight bracket from 80,000 Ib. to 120,000 Ib. and a temperature bracket 
32) from I.S.A. to 1.$.A.+30°C., then the reference climb weight and temperature 


should be 100,000 Ib. and I.S.A.+ 15°C. respectively. 
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The calculations should then be made for (X —D)V/W for the reference climb 
weight and temperature selected as just mentioned, and plotted against energy 
height. Such a plot may be linear or very close to linear up to an energy height at 
which (X—D)V/W is about 300 or 400 ft./min., in which case, for practical 
purposes, a straight line may be drawn through these points and the relation may 
be expressed as 


where B and C are constants. 


The time integral then becomes 
dH.» 
t= | ‘ . (36) 


As in Part I the concept of the relative energy height is used. Integral (36) can 
be evaluated in the usual way. 


When a plot of (X—D)V/W against energy height for a reference climb 
weight and temperature is not linear, as assumed in equation (35), the relationship 
could be expressed as 


where a, b and c are constants. The constants a and b can be determined by 
differentiating equation (37) with respect to energy height. 


The time integral becomes 


dA. 
3 


0 


Integral (38) may be evaluated in the same manner as equation (12) in Part I. 


The distance and fuel integrals in both these cases may then be evaluated 
either mathematically, as has been done in Part I if fuel flow and climb speed can 
be expressed approximately as linear functions of energy height; otherwise 
graphically, as in the conventional method. 


Once the constants a, b, c or B and C are determined for a reference climb 
weight and temperature then the effect of a change of climb weight and/or 
temperature from the reference case on these constants can be evaluated approxi- 
mately by the formulae given in Appendix I. Thus, once a relationship between 
(X —D)V/W and energy height is established in the form of either equation (35) 
or (37) for a reference climb weight and temperature, a similar relationship for 
other climb weights and/or temperature conditions can be established in the same 
form without any further calculations for (X —D)V/W for those particular cases. 
The time, distance and fuel integrals can be evaluated as already discussed. 
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8. Discussion 


In equations (6) to (10) the constants n, g, h, m and r which are slopes of 
thrust, fuel flow, climb speed (T.A.S.), drag and climb weight curves with energy 
height are assumed to be independent of temperature. Similarly the constants 
nw, q’, h’, m’ and r’ which are slopes of the thrust, fuel flow, climb speed, drag 
and climb weight curves with temperature are assumed to be independent of 
energy height. These criteria will not rigorously be met in practice. However, if 
the variations of these slopes are not appreciable, a mean value may be assumed 
in the calculations. 


For better accuracy, however, a slightly different procedure described under 
Case I (A) may be followed. That is to say, all the five equations from (17) to (21) 
should be established for different climb weights and temperature conditions to be 
considered. Then by using the appropriate formulae the time, distance and fuel 
integrals should be evaluated. 


Where there is a change of slope of speed or thrust or fuel flow curve with 
energy height the method could be applied in two stages. For example, it may 
be desirable for a turbo-jet aircraft to make a climb at certain speed and engine 
I.p.m. up to a certain altitude and then change the climb speed and/or engine 
setting. This change in climb technique at certain altitude, demanded by certain 
reasons, will cause changes in the slopes of the curves of speed, thrust and fuel 
flow against energy height above that particular altitude. In these circumstances it 
is desirable to evaluate the time, distance and fuel integrals in two stages. 


The accuracy of the results by this energy method depends entirely on the 
degree with which the basic assumption of linear variation of thrust, fuel flow, 
speed and drag with energy height is met. It is therefore mandatory that the 
variation of thrust, fuel flow, climb speed and drag (based on constant weight for 
this purpose) with energy height should be examined before a calculation is under- 
taken. The degree of deviation of these quantities from the linear assumption will 
give an approximate idea of discrepancy between the estimated results and those 
which may be expected from flight tests. 


Good approximate results with some saving in time required for climb calcula- 
tions may be obtained by using a constant climb weight (an average weight 
during climb). 


Moderately satisfactory results may be obtained by using the method described 
in Part II. This method is simpler and quicker inasmuch as it is necessary to 
establish a relationship of (¥-—D)V/W against energy height only for a reference 
climb weight and temperature, in the form of either equation (35) or (37). Then 
that form of equation can be extended to other weights and/or temperature 
conditions without actually doing the calculations for and making a plot of 
(X-D)V/W against energy height for those particular conditions. Once again 
the degree of accuracy of the results will depend on the degree with which the 
assumptions made in deriving equations (47), (49), (50), (51), and so on, in 
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Appendix I are met. By virtue of the approximations and assumptions made in 
deriving equation (35) or equations (47), (49), and so on, just mentioned, a high 
degree of accuracy of the results should not be expected. Nevertheless this method 
will give satisfactory results with considerable saving in the time required for 
climb calculations. 


9. Conclusions 


As a result of the investigations carried out in developing and applying the 
energy approach method described in this paper to the climb performance 
estimation of a civil turbo-jet transport aircraft, the following conclusions are drawn. 


(i) The application of the energy approach method to the climb performance 
estimation of a turbo-jet transport aircraft for various climb weights and tempera- 
ture conditions will prove satisfactory from the points of view both of the accuracy 
of the results and of saving some of the considerable time needed for such 
calculations by the conventional method. 


Some idea of the time saved by this method during the calculations of a 
group of three different climb weights and temperature conditions can be obtained 
from the following table, which is based on the author’s experience. 


Conventional method Energy approach method 
Part I Part II 
2 : 1 


However, if the climb calculation is desired for only one climb weight and 
temperature condition, then there is very little to choose between the energy 
approach method and the conventional method. The conventional method in this 
case is quite satisfactory. 


(ii) The energy approach method will prove useful both in the preliminary 
design stages and in the development stages of a turbo-jet transport aircraft as a 
quick means of studying the effects of small changes in one or more of the controlling 
parameters of climb performance, such as thrust, drag characteristics, and so on, 
on the climb performance of an aircraft under various climb conditions. 


(iii) The method developed in Part II is simpler and quicker than that in 
Part I; as a result its application will save still more of the time required for such a 
calculation. 


(iv) The operational and development sections of the airlines should also 
benefit by this method. 


(v) It is thought that the application of this method to the descent performance 
estimation, as well as to the calculation of the drift-down after an engine failure, 
will also prove satisfactory. 
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Appendix I 


EFFECT ON THE SLOPE OF THE CURVE OF CLIMB WEIGHT AGAINST ENERGY HEIGHT 
OF A CHANGE IN (a) CLIMB WEIGHT AND (b) TEMPERATURE 


(a) At I.S.A. conditions (or at the reference temperature considered, that is at 
AT=0) the relation between the weight in climb and energy height can be expressed 
approximately as 


The decrease in weight during climb is equal to the fuel consumed in climb. 


dQ= -dw= ~FxAt 


__ 
(X¥-D)V/W 


dw F 


where F=mean fuel flow 
At=time to climb energy height, dH... 
(X — D) V/W=mean rate of climb during height, dH,,. 


By equation (41), r is a function of fuel flow and rate of climb, which in turn are 
functions of energy height. For a first approximation an average value of r may be 
evaluated at 67 per cent of the expected cruise height at a reference climb weight. 


The drag during climb, which depends on the weight, may be evaluated by successive 
approximations at the same height at which r is to be evaluated. In practice two 
successive approximations are found to be satisfactory. As thrust, fuel flow and climb 
speed, by assumption, are independent of climb weight, r could be evaluated by using 
equation (41) for different climb weights. 
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So that equation (41) may be written as 


le - D) V/ W 67 per cent cruise height ( ) 


That is, F, X, D and V are to be evaluated at 67 per cent of the expected cruise 
height at a reference climb weight. 


(b) Temperature effect, like weight effect on r, may be evaluated in a similar 
manner, so that 


43 
le a D) V/ W 67 per cent cruise height ¢ ) 


where F, X, D and V are evaluated under appropriate temperature conditions at 67 per 
cent of the expected cruise height. 


2. EFFECT OF A CHANGE IN CLIMB WEIGHT ON CONSTANTS 4a, b, c, B AND C IN EQUATIONS 
(35) AND (37) aT I.S.A. CONDITIONS (OR AT REFERENCE TEMPERATURE 
CONSIDERED) 


For the purpose of calculating the effects of change in climb weight and temperature 
on the constants a, b, c, B and C in equations (35) and (37) the following simplifications 
and approximations are made :— 


(i) The thrust and climb speed (T.A.S.) are approximately linear with energy 
height. (If however thrust and climb speed are expressed as some other 
convenient function of energy height, a similar treatment on the lines 
mentioned in what follows can easily be developed for those cases). 


(ii) The weight of the aircraft during climb is held constant. 


By virtue of these assumptions and approximations, equations (17) to (21) in the 
main text can be used. The basic expression for (X - D) V/W then becomes 
1 


V 
(X-D) = wl 


V, (X,- D,)+ Her [V, (n- m) +h (X, D,)] + hHer? (n- m)}. (44) 


Constants C in equation (35) and c in equation (37) are given by 
and B in equation (35) is given by 


d 


B= 
dH.» 


[x ~D) = [V, (n- (46) 


Assuming m to be constant with weight and 


Ww? 
b,2ne4pV? 
=G+RW?, 
1 
where G=Cy,+pV Bes R= by? pV?” 
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equation (46) becomes 
B= (n- m)+h(X,- G- RW?’)}. ‘ . (46a) 


The effect of change of weight on the constant B can be evaluated by differentiating 
equation (46a) with respect to weight. 


d B 
Equation (44) may be rearranged as follows: — 


V 
(X - D) W bH.p+c 


Equation (48) is the same as equation (37). 


The effect of change of weight on constants a, b and c in equation (48) may be 
evaluated by differentiating these constants with respect to weight. 


d a 
Thus aw (a)= - W . (49) 
d b 
aw (b)= - (2nr +. (50) 


Equations (47), (49), (50) and (51) are obviously valid for assessing the effect of only 
small changes in climb weight of the aircraft. For large variation in climb weight, 
however, the modified formulae given in what follows should be used. 


Let W, be a reference weight and W,—W, + AW, where AW is 5 per cent or more 
of the reference weight. From equation (46a) 


B W 
B, - B,=dB= - (W, W.) {= +Rh(1+ : . (52) 
1 2 1 2 W, W, 
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A similar approach will modify equations (49), (50) and (51) to 


a, —a,—da= —(W,—W,) 


Ww 
c,—c,=de= —(W, {© Rv, (14 
2 2 


3. EFFECT OF CHANGE IN TEMPERATURE ON CONSTANTS a, b, c, B AND C IN 
EquaTIons (35) AND (37) AT CONSTANT WEIGHT 


By virtue of the assumptions and approximations made in the previous section, 
equations (6) to (9) in the main text can be used. Thus the basic expression for 
(X - D) V/W becomes 


(X —D) = { (x, +n WAT) + Hep [ (n—m) HAT) + 


Differentiating equation (56) with respect to AH,,, the constant B in equation (35) is 
obtained as 


d V 1 
+h(X,—D,)+AT (n—m)+h(n'— (57) 


As AT is a relative temperature with respect to the I.S.A. conditions or the reference 
temperature considered, the difference in the value of B, (say) at the I.S.A. condition (or 
at the reference temperature) and B, (say) at a temperature which is equal to I.S.A. 
(or the reference temperature) +A7°C is obviously given from equation (57) as 


1 
B, —B,=dB= +AT < 8) 
At the I.S.A. conditions (or at the reference temperature), that is AT=0, dB=0. 
Now equation (56) could be rearranged as 
V 
D) W = + bH +C. 


h(n—m) 
where a= Wr 


b= {v, (n—m)+h(X,—D,)+AT (n—m)-+ h(n — my} 


c= {(x,-D,) Vit AT LA — m’) +h’ (X,—D,)]+ AT?h’ —m)} 
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Equation (59) is the same as equation (37). 


The difference between a,, b, and c, (say) at the I.S.A. condition (or at the 
reference temperature) and a,, b, and c, (say) at a temperature which is equal to the 
1S.A. (or the reference temperature) +A7°C is obtained from equation (59) 
as follows :— 


a,—a,=da=0. 


b,—b,=db=+4AT { [h’ (n—m)+h my} 


C,—C,=dc= + { AT [V,, (n’ — +h’ (X, — D,)] + (n’— m’) 


If the constants h, n and m are not independent of temperature, the constants a, b, c, B 
and C ought to be determined for each temperature condition to be considered. If, 
however, h, n and m are very nearly independent of temperature (this will be the case 
for a small variation of temperature from the reference case), then equations (58), (60), 
(61) and (62) may be used to evaluate an effect of change in temperature on a, J, c, 
B and C. 


m’ in equations (58), (61) and (62) is usually very small in comparison with n’ and 
hence may be neglected. 


Appendix II 
ILLUSTRATIVE EXAMPLE 


For a turbo-jet transport it is required to estimate the time needed, distance covered 
and fuel consumed in a climb from 5,000 ft. to 35,000 ft. pressure altitude under the 
following conditions : — 


(i) climb speed: constant E.A.S.; 260 m.p.h. 
(ii) engine speed: maximum continuous r.p.m. 
(iii) climb weight: 80,000 Ib., 100,000 Ib. and 120,000 Ib. at LS.A.+ 10°C. 
and (iv) at temperature I.S.A.-—5°C. and I.S.A.+25°C. for 100,000 Ib. 


The reference climb weight and temperature in this case are 100,000 lb. and 
LS.A.+10°C. respectively. It is thus intended to cover a climb weight bracket of 
20 per cent of the reference weight and the temperature bracket is equal to the 
reference temperature + 15°C. 


From the calculations and plots of thrust, fuel flow and climb speed (T.A.S.) against 
energy height under the temperature conditions just mentioned, the constants in 
equations (17), (18) and (19) are determined by assuming a linear relationship of those 
quantities with energy height. (See Table I.) 


Using equations (42) and (43), first the average values of the slope r of the curve of 
climb weight against energy height were determined and then the average values of 
the slope m of the curve of drag against energy height were established. The results 
of the calculations are given in Table II. 
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TABLE I 
CONSTANTS IN EQUATIONS (17), (18), (19) 
; 
(ft. ft. ) (Ib. (Ib. Ib. /sec. 
af | ft. ) 
1.S8.A.—5°C | 396 8-10 x 10- 3 | 17: 8 x 10° 255 4°550 —7°41 x 10-5 
1.S.A. + 10°C. 408 7°89 x 10-3 x 108 —0-200 4:085 — 6°40 x 10-5 
1S.A. +25°C. | 418 7:55 x 10-3 | 13-9 x 10 —0:°173 3-700 — 5°75 x 10-5 
V,, X, and F, correspond to 5,000 ft. (H, ) from which the climb is commenced. 
TABLE II 
VALUES OF r IN W=W,+rH.,; AND m IN D=D,+mH,, 
Climb weight: Ib. 80, “100, 000 120,000 
r ft. —0:063 — 0-087 —0°121 
m: |b./ft. 54 x 3 —2:°65 x 10-3 —4:4x 10-3 
100,000 Ib. 
“ag S. A. Ss. A.+ +10°C. 1.8.A.+25°C. 
ib:/ft. —0°8 — 0-087 —0°102 
m: Ib./ft. —2°44 x 10-3 —2°65 x 10-3 —3-1x 10-3 


Part I. Having thus established all the constants in equations from (17) to (21), the 
time, distance and fuel integrals were evaluated by using equations (23), (24) and (25) 
respectively. The results of the calculations are presented in Table IV. 


Part II, The actual plot of (¥- D)V/W against energy height at the reference 
weight and temperature (100,000 Ib. and I.S.A.+ 10°C.) was not quite linear, but an 
approximate linear relationship was assumed so that 


=C + 


where C=42:7 ft./sec. 
B= — 6-47 x 10-‘ ft./sec./ft. 
Constant C here is not the actual value of (X-D)V/W when H.,=0, but an 


intercept on the ordinate of the curve of (X¥—D)V/W against H,,, assumed a 
straight line. 


The effect of the change in weight from 100,000 Ib. to 80,000 Ib. and 120,000 Ib. on 
the values of the constants B and C was calculated by using equations (52) and (55). 


The temperature effect on the value of B was evaluated by making a straight- 
forward calculation using equation (57) at a mean energy height, as equation (58) could 
not be used because A, n and m are not independent of temperature (see Table I). 


As C in this case is an intercept on the ordinate of the curve of (X - D) V/W against 
energy height, it was corrected for different temperature by using equation (62). 


The values of B and C are collected in Table III. 
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TABLE III 
VALUES OF B AND C IN(X—D)V/W=C+B H,, 


1.S.A.+10°C. 

"Climb weight: Ib. 80,000 100,000 420,000 
B: ~7:51x 10-4 ~6-47x 10-4 ~5-85 x 10-4 
C: ft./sec. 56°2 42:7 33-2 

Temperature 1.S.A.—5°C. 1.S.A.+ 10°C. 1.8.A.+25°C. 
B: —5-96 x 10-4 ~6-47 x 10-4 — 4-96 x 10-4 
ft. 49-4 42:7 35°5 


Using the values of B and C given in Table III the time, distance and fuel 
integrals were evaluated and the final results of the calculations are presented in 


Table IV. 


For convenience the final results of the calculations for time are given in minutes 


instead of seconds and for distance in statute miles instead of feet. 


SUMMARY OF THE RESULTS 


The results of the climb calculations in terms of the time needed, distance covered, 
fuel consumed and air miles per pound of fuel in climb from 5,000 ft. to 35,000 ft. 


pressure altitude, in the illustrative example, are summarised in Table IV. 


TABLE IV 
1.S.A.+ 10°C. 

80,000 100,000 120000 

By By By By By By 
method method method method method method 
in in in in in in 
Part I Part Il Part I Part Il Part I Part Il 
Time: min. 14-75 14-75 20°6 20°6 28:8 28°8 
Distance : 
statute miles 96 95 133 134 189 187 
Fuel: Ib. 2,480 2,470 3,430 3,424 4,790 4,738 
Air miles per Ib. 
of fuel x 102 3°87 3°85 | 3°88 3-91 3:95 3-95 
100,000 1b. 

Temperature 1.S8.A.—5°C. 1.S.A.+ 10°C. 1.S.A.+25°C. 
Time: min. 15-2 14:35 20°6 20-6 29-75 26°5 
Distance : 

statute miles 95 89-4 133 134 195-5 178°5 
Fuel: Ib. 2,872 2,825 3,430 3,424 4,220 3,935 

Air miles per Ib. 
of fuel x 102 3°31 3°17 3-88 3-91 4-62 4°54 
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COMMENTS ON THE RESULTS 


The method described in Part II of the paper, in spite of its approximate nature, 
gives results in this case which are satisfactory when compared with the results obtained 
by the method described in Part I. The maximum discrepancy between the results of 
the two methods occurs for a weight of 100,000 Ib. with I.S.A. + 25°C, where it is about 
8-7 per cent in distance covered and about 6-7 per cent in fuel used during climb. It 
is interesting to note that the discrepancy in the results obtained by the two methods, 
whatever its magnitude may be, is of the same sign both for the distance covered and 
fuel consumed during climb. For this reason the air miles per pound of fuel during 
climb are about the same by both the methods. The discrepancy in the air miles per 
pound of fuel during climb in the case given (I.S.A.+25°C) is only 1-7 per cent, the 
maximum being 4-2 per cent at I.S.A.—5°C. 


The comparison which was made separately of the results of calculations by these 
two methods with those obtained by the conventional method has also shown the same 
tendency for fuel and distance discrepancy of the same sign. 
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The Natural Frequencies of Thin Skew Plates 


R. K. KAUL and V. CADAMBE 


(National Physical Laboratory of India) 


Summary: Using Rayleigh’s principle, the natural frequencies of thin isotropic 
rhombic plates, with three possible combinations of boundary conditions 
obtained by combining clamped-clamped and clamped-supported edge con- 
ditions, are determined in Part I of this paper. To introduce constant limits 
of integration, non-orthogonal co-ordinate systems are used and the wave 
shape for the vibrating plate is approximated by using normal functions 
representing mode shapes of corresponding bars. To estimate the accuracy 
of these eigenvalues, Kato’s theorem is used and the lower bounds for the 
natural frequencies are determined in Part II of the paper. It is also shown 
that normal beam functions are not generally suitable for the determination of 
eigenfrequencies of skew plates with large skew angles. 


1. Introduction 


The study of natural inextensional vibrations of thin oblique plates, with 
particular reference to the determination of eigenfrequencies and mode shapes, is a 
problem of considerable interest. Other methods are available but Rayleigh’s 
minimal principle”? finds great favour in spite of certain limitations. One of the 
difficulties is the suitable selection of a minimising sequence and therefore quite 
often, to get a closer approximation to the natural frequencies, it it necessary to 
resort to the “ device of variable parameters ” apparently due to Rayleigh, although 
sometimes attributed to Ritz’. 


Rayleigh’s energy principles shows that first order errors in the tentative mode 
of vibration produce only second order effects in the estimate of natural frequency 
and, since it is not possible to admit all possible functions, the Rayleigh’s ratio is 
therefore no longer the exact solution A,* but only a close upper bound to the 
natural frequency. In a number of cases this principle leads to surprisingly close 
results; but, in cases where exact values of the frequencies cannot be calculated, a 
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criterion is needed to judge the accuracy of the results obtained. The only justifica- 
tion is then to calculate the close lower bounds to the natural frequencies, a know- 
ledge of which would be as useful, for practical purposes, as an exact value of the 
frequency. 


By using a condition of “ weak minimum ”, Weinstein has shown how lower 
bounds to a clamped plate problem can be obtained. But recently Kato in a 
singularly simple and elegant theorem gave a method, which applies not merely to 
fundamental mode, but to any overtone of a vibrating system with any number of 
degrees of freedom, finite or infinite. Temple‘) (1952) extended Kato’s theorem to 
vibrating systems with both elastic and inertial coupling and further obtained 
generalisations of his results. 


In Part I of this paper, using Rayleigh’s principle, approximate frequency 
equations are derived and eigenvalues determined for three different boundary 
conditions of skew plates when vibrating in modes of the type m/n and (m/n+n/m). 
In Part II, Kato’s theorem is used for determining the close lower bounds for which 
the upper bounds are provided by Rayleigh’s principle. 


NOTATION 
reduced elastic strain energy of the bending of the plate 


reduced kinetic energy of the plate vibrating harmonically 


Young’s modulus and Poisson’s ratio of the plate material 
respectively 


flexural stiffness, Eh*/ {12 (1 —v?)} 


”frestangua and oblique co-ordinates, respectively 


+ 


7, constant mass per unit area of the plate of uniform thickness h 
® natural circular frequency, 27f 

A,? exact eigenfrequency, (c,/D,) cos*a 

Pr” Rayleigh’s ratio 


®,,(), ¥,.(Q) normed-orthogonal bar eigenfunctions representing normal modes 
of vibration of a uniform bar 


=0, for 
Kronecker’s delta 

=1,form=n 
measure of accuracy 
residual inertial force 
symbols used in Kato’s lemma 


two eigenvalues defined in Part II 
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VIBRATION OF THIN SKEW PLATES 


PART I 


ON THE CALCULATION OF UPPER BOUNDS TO NATURAL 
FREQUENCIES 


2. Basic Theory 


The equation of motion of a thin isotropic plate is governed by the fourth order 
partial differential equation 


=0 in domain D (1) 


where w, is the displacement at any point (x, y) at time ¢. It is assumed that the 
motion is simple harmonic and therefore 


w, =w sin (wf + 4), 


where w satisfies the boundary conditions of the plate. The formal approach to the 
problem consists in selecting for w a double infinite series, each term of which 
satisfies equation (1) and some of the boundary conditions. Then, by judiciously 
combining an infinite number of these, the remaining boundary conditions can be 
satisfied. This procedure, first adopted by Taylor’ in the analysis of statical 
problems of buckling of a clamped rectangular plate, was later used by Tomotika‘”’ 
for a similar problem in vibrations of plates. The procedure leads to an infinite 
determinant, the roots of which give the natural frequencies. Normally, by taking 
a few rows and columns, the roots converge, but the process is fairly laborious. 


Using Rayleigh’s minimal principle the natural frequencies of a vibrating system 
can be quite closely estimated. If the reduced potential and kinetic energies of a 
continuous system be represented by V,, and T,, then Rayleigh’s principle states 
that for a tentative mode w, pr?=V,,/T, is stationary, i.e. first order errors in the 
tentative mode produce only second order effects in the frequency. When no exact 
solution can be found, w can be approximated in the form 


+ %,W, + 


in which w, satisfies the artificial boundary constraints of the problem and the 
constant co-ordinates «, are to be so adjusted as to render V,,/7, a minimum. 
Differentiating with respect to these arbitrary parameters and equating to zero, the 
general frequency equation is 


0a, {Vin pr Tw} =0 
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For thin plates the quadratic functionals V,, and T,, are defined by the following 
bi-linear expressions : — 


| - (se) } ] dx dy 


(3) 
and | dx dy 


For plates of polygonal shapes, when w vanishes on the boundary, it can be 


shown that 
02w } 
ox dx dy= 


and therefore V.= dx dy. 


Hence the eigenfrequencies for such plates are given by the simpler equations 


Om, 


-pp?w*}dxdy=0 . ; ‘ (4) 


Fig. 1. 


Thin plate of oblique plan form. 


For a plate of oblique plan form, as shown in Fig. 1, the co-ordinate transformation 
relations are 


€=x-ytana 


(=v sec 2, 
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and therefore = ay =SeC 2 — sin 


and Jacobian J y ) =COS 


Therefore the equivalent frequency equations in terms of the skew co-ordinate 
system (€, ¢) are 


pxiw* | dé 


where pr? =o? — cost a, 
D, 


This leads to a system of n simultaneous linear homogeneous equations in n coeffi- 
cients and the eigenvalues p,z,”. For a non-trivial solution to exist it is necessary 
that the determinant of the system of equations be zero; the successive roots in the 
ascending order pr,” < pro” S pr;? S »<Prn” are the discrete eigenvalues of 
the continuous system. Since in the selection of an assumed form for a corresponding 
mode of vibration all possible modes are not admitted, the estimate of eigenvalue 
on the basis of the stationary property will always err on the side of excess, and 
therefore Rayleigh’s ratio pp? =>A*. The order by which Rayleigh’s estimate 
approximates the exact frequency depends however on the selection of the plate 
deflection function and the closeness with which it represents the actual wave form. 
Though it is possible to improve the approximation and refine the method into a 
rapidly convergent procedure, the detail of the numerical computation increases 
inordinately with each improvement of approximation and therefore becomes more 
and more laborious. 


3. Bar Eigenfunctions 


For the skew plate with edges parallel to € and (¢ co-ordinates, the series 
approximation for the plate deflection w is taken in the form 


where 2,,, are the arbitrary parameters consistent with the required normalisation 
condition, and ®,,(€) and W,(¢) are sets of “ admissible ” functions. These functions 
admissible for competition must explicitly satisfy the “ artificial boundary condi- 
tions”, since fixed boundary values of these functions, or their first derivatives which 
imply slope, correspond to rigid constraints on the boundary of the plate. On the 
other hand, natural boundary conditions need not be stipulated in advance, as these 
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express the equilibrium of the plate along the boundary, and are implied in the 
variational problem’. In the present case a very good approximation to the 
admissible set of deflection functions is the choice of eigenfunctions defining normal 
modes of vibrations of corresponding bars and which satisfy the required boundary 
conditions. These normed-orthogonal bar eigenfunctions are: — 


(i) Clamped-clamped bar 
sin {km(E-@/2)} sinh {kn (€ - a/2)} 

1 cos {km(€-a/2)} cosh {k,, (€- a/2)} mr 


sin? —— 
Va cos (k,a/2) cosh (k,a/2) 2 


+ 


for m=1, 2, 3, ‘ (6) 
where k,,a is the m™ positive root of the transcendental equation 
tan (k,a/2)=(- 1)” tanh (k,,a/2) 
or its equivalent cosh kna=sec kna. 


(ii) Clamped-supported bar 


sin {km(€-a)} sinh {kn - a)} } 


cos cosh Kna 


for m=1, 2, 3, 
where k,,a is the m™ root of the transcendental equation 


tan k,,a=tanh k,,a. 


The functions are not only continuously differentiable in the range 0 <= <a 
but also satisfy the homogeneous differential equation 


and the terminal conditions of corresponding uniform bars vibrating harmonically. 
Further the set of functions selected are normed-orthogonal within the range of 
integration and satisfy the condition 


#, ©, © 


where 54,, is the Kronecker’s delta. 
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The second derivatives of these functions are also orthogonal and satisfy 
the relations 


0 


where k,a is the r bar eigenvalue. These ortho-normal bar eigenfunctions, in 
addition to the usual properties of Fourier series also possess the property 


a a 


0 


Similar properties also hold good for the function W,(¢). 


4. Approximate Frequency Equation (Rayleigh’s) 


Using Rayleigh’s principle, the tentative shape for the type mode m/n is taken 
in the simplest form 


(€) ‘ ‘ ‘ (9) 


where ®,,(€) and W,(¢) are the m™ and n‘ eigenfunctions of corresponding bars in 
the €- and (-directions. On substituting for w in the frequency equation (5) and 
using the properties of the normed-orthogonal functions described earlier, the follow- 
ing equation for the Rayleigh’s ratio is obtained 


> ab? \ (k,,a)* + (k,,b) + 


+ 2a°b’ (1 +2 sin? 2) | Vn (E) Pm” (E) Vn Un” (Q at} (0) 
0 


Besides modes of the type m/n, i.e. m half-waves along the €- direction and n 
half-waves along the ¢- direction, the nodal patterns of which are lines parallel 
to the co-ordinate directions, there also exist normal modes of the type (m/n+n/m) 
when the frequencies of the modes m/n and n/m are identical. The nodal patterns 
for these are formed from the simultaneous existence of two component vibrations 
m/n and n/m at a point and the resultant amplitude is governed by their in-phase 
or out-of-phase relationship. The two new modes (m/n+n/m) have two discrete 
eigenfrequencies and the nodal patterns are also discrete. If the wave form for 
such a mode is represented by 


November 1956 343 


he 
al 
ry 

) 
ab 
) 


R. K. KAUL AND V. CADAMBE 


then the eigenfrequency p,” for the distinctly discrete set (m/n+n/m) is given by 


> 


1 
Pa = [a {(pma)* +(pna)*} + {(Gmb)* + (qub)*} + 


+ 2a?b? (1+2 sin? {Pn () Pn” © O+ 
) 


+P, (€) Pr” (€) Vn Un” Q42Pm Pn” (E) Un Vn” O} | (12) 


Expressions for the various definite integrals which occur in the frequency equations 
are given in the Appendix. 


(a) Oblique plate clamped on all edges 


For clamped-clamped skew plates, the same bar eigenfunction (6) is used along 
the two co-ordinate directions. In the frequency equation (12) the integral 


6 


0 


is zero when (m+n) is odd, and therefore in this case discrete modes of the type 
(m/n+n/m) exist only for those combinations when (m-n)= +2, +4, '+6, 

For other combinations of m and n, only modes of the type m/n exist, the eigen- 
frequencies of which are given by equation (10). For a square plate or when the 
skew angle is small (2 < 15°), the nodal patterns corresponding to modes m/n are 
straight lines parallel to the edges of the plate. But for modes of the type 
(m/n+n/m) it is not so and the amplitude of any point, depending upon the phase 
relationship, is the sum or difference of the amplitudes of that point for the 
component vibrations m/n and n/m. 


(b) Oblique plate clamped on three edges and supported on the fourth 


Both eigenfunctions (6) and (7) are necessary in this case to represent the plate 
deflection wave form. The function ®,,(€) represents the deflection of the plate 
in the clamped-clamped direction and the function ¥, (¢) in the clamped-supported 
direction. Although in the frequency equation (12) the integral Imn-Jmn has definite 
positive values for those combinations of m and n when (m-n)= +2, +4, 
associated modes of the type of (m/n+n/m) cannot exist, since for this unsym- 
metrical set of boundary conditions the eigenfrequencies of modes m/n and n/m 
are different. Therefore for all combinations of m and n, only modes of the type 
m/n will exist. 


(c) Oblique plate with two adjacent edges clamped and the other two freely 
supported 


The clamped-supported bar eigenfunction (7) is used to represent the deflection 
of the plate in this case. The integral mn Or Jmn has a definite value for every 
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combination of m and n, (mn). Further the boundary conditions are symmetrical 
along the diagonal of the plate and therefore eigenfrequencies of modes m/n and 
n/m will be identical, in which case compounding of modes will take place, leading 
to discrete modes (m/n+n/m). In this particular case the eigenfrequencies of the 
two associated modes (m/n+n/m) are very nearly equal. This is because the 
contribution of the integral /,,,-J;,., in the frequency equation (12) is very small. 
Thus, if the mode is wrongly assumed to be pure m/n and the effect of n/m is 
neglected, the error introduced in the results will be negligible. 


For each of the three sets of boundary conditions just discussed, Rayleigh’s 
ratio is calculated for both modes of the type m/n and (m/n+n/m). Plates having 
skew angles of 0°, 15°, 30° and 45° and lateral dimensions (a x a) are considered 


and the successive eigenvalues in the ascending order 
— 


Pri BS Pre S Pr3 S 


are given in Table I. 
PART II 


ON THE CALCULATION OF LOWER BOUNDS TO NATURAL 
FREQUENCIES 


5. Basic Theory 


From a practical standpoint (when an amount of uncertainty is involved in 
the physical data) it is essential to know the margin by which the frequencies err 
on the side of over-estimation. Rayleigh’s principle therefore suffers from a serious 
drawback in that it provides only an upper estimate for the frequencies. Rayleigh 
himself has given two “ Comparison Theorems ” which provide the lower estimate, 
but with far less accuracy than his energy method provides the upper bounds. The 
singularly simple theorem of Kato answers the requirement in that it provides really 
close lower bounds to the natural frequencies. The original results of Kato are 
restricted to systems with no inertial coupling and do not apply to general 
Lagrangian systems. These were recently extended by Temple’ to systems with 
both elastic and inertial coupling and are used here in the determination of lower 
bounds for the natural frequencies of skew plates for which the upper estimate (in 
Part I) is based on Rayleigh’s principle. The theory of the principle is taken mutatis 
mutandis from Ref. 4 and therefore only portions relevant to the problem under 
consideration are discussed here. 


Let 2” and 6? be two eigenvalues corresponding to two distinct normal modes 
of vibrations (which by convention are normed), and selected such that 


while no eigenvalue 4,” of the system lies in the interval 
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To maintain oscillations in the trial mode w for which the approximate eigenvalue 
iS pr”, a residual frictionless force R will come into play and therefore, for the 
trial mode w, 


Similarly, if the residual forces corresponding to two eigenvalues <* and f? be 
A and B, then 


(A. B)=(V'w - B?w) 
= 2*) w} {R+ (pr? - 8’) w} 
and since external forces are supposed to do no work on the system, we have 


(A e By=R ~R+ (pr? 2°) = B?) w.w a (14) 


Further, in accordance with Kato’s lemma, by using the theory of normal modes 
and principal co-ordinates, it can be shown that 


K,.=(A. B)= (A? - 27) (A,2 - B?). ; (5) 
But 2? and £? have been chosen so that no eigenvalue 4,” lies in the interval 
Therefore (A,? — z*) and (A,? - 8*) will always have the same sign, making 
K,=(A.B)>0. 


Hence, a fortiori, from equation (14), 
(pr? — 2°) (B? Pre) S 


where ¢,7=(R.R)/(w.w), 


and in rectangular co-ordinates 


In applying this principle for determining the lower bounds to an eigenvalue 
(say), for which the close upper bound is >A,*, we take z?=A,? and =p? 
where p’ is the smallest eigenvalue greater than A,* and a lower estimate to A,’. 
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Although such a choice implies multiplicity of A,?, at the same time it fulfils the 
basic condition of Kato’s lemma that no eigenvalue A,” lies in the interval 
a? <A,?< Since pr,? >A,?, we have 


< 
{ (8? — pr’) } 


and therefore the bounds to A,” are given by the “ bracket ” 


{ px? - <p,’ 


6. Application to Skew Plates 
In the skew co-ordinate system, the equation of motion of a thin plate, 


corresponding to equation (1), is 


4 sin }+4sintz -atw=0 . (18) 


where ( a and =? cos* a, 
> 0 


For Rayleigh’s ratio pz”, determined on the basis of an assumed mode shape, the 
measure of accuracy is 


ab 
+4 


For the type mode m/n, we take the simplest approximation (equation (9)) for 
the mode shape and, since ®?,,(€) and W,(¢) are normed-orthogonal functions, we 
have, on substituting in equation (19), 


= (Km* + kn*) (Kint + kn* + prt +4{(1 + 2 sin? 2)? +8 sin? a} k,,*k,* + 


4 
b 

+4(1+2 sin? 2) (k,,4+k,* pr) Vn(E) (Q) 
0 


b 

16 sin? 
0 


(20) 
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where, using the terminal conditions, the transformations are 


a a 


and 4%) dé = [ kn! Dy(E) dE 
0 0 0 


for clamped-clamped and clamped-supported functions. 


Equation (20) holds good for the general case of oblique plates with skew 
angle 2, vibrating in mode m/n. When the plate vibrates in an associated mode 
(m/n+n/m), the simplest approximation for the mode shape is taken in the form 


Substituting for w in equation (19), the measure of accuracy for the particular case 
of a square plate is 


= (Kin! + (Kent + kn! + + + 


P,, (x) ?,,” (x) dx, (x) ®,” (x) dx 


Jan (y) (y) dy, Pn (y) (y) dy. 


Corresponding expressions for <,? can similarly be worked out when the trial 
mode is made up of a number of vector components with arbitrary parameters ~,.. 
For determining the lower bounds, it is now only necessary to know §? (a lower 
estimate to the next higher eigenvalue A,”). This can be calculated by using 
Rayleigh’s Comparison Theorem’, and for a continuous system typified by 


Tw-Vw 


(where V and T denote the bounded operators which specify the potential and 
kinetic energies) it is given by the smallest of the quotients {Vw/Tw}, for a 
necessarily positive product (w. Tw). As an alternative, Rayleigh’s ratio p,’” for 
the next eigenvalue may be used as an approximation. Kato’s inequality (17) now 
yields the lower bounds. 


7. Comparative Results and Discussion 


On the basis of the theory developed in the previous section, lower bounds to 
the eigenvalues of a skew plate are determined, for which the upper estimate on the 
basis of the Rayleigh’s principle is obtained in Part I. Taking simplest tentative 
modes, the measure of accuracy is calculated in each case from equations (20) or 
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‘TABLE Il 
LIMITING BOUNDS OF RHOMBIC SKEW PLATES* 


k= { , Kato’s lower bound. 
Py?) 


cos?z, Rayleigh’s upper bound. 
0 


Percentage 
Mode | Lower and upper bounds | Mean 


deviation 
type Px S*, SPr from mean 


(a) Chemport clamped skew plate 


3533322 3610868 35-72095 
71°76754 73-73720 72:75237 
104-98783 108-84992 106-91886 


3469011 36°66593 35°67802 
63°68607 74:°75947 69:22277 


32°95941 38-14697 39°55319 


30°63837 40: 08173 35:36005 


on three and on the fourth 


3146043 31:95364 |  31-70704 
62:22732 63-65944 62:94338 


31°46798 32:°54105 | 32:00451 
60°88146 64°76057 62°82102 


30-35069 34:09421 | 32-22245 


| 29-46388 | 32-78597 


(c) € lamped 0 on two adjacent edges alt supported on the other two 


1/1 2622513 27 19478 26°70995 
| 5940657 60°69028 60°04842 


| 
| 
| 
| 
| 
| 


1/1 | 24:91261 27:83775 26°37517 


1/1 21°45018 29-52310 25°48664 


* Common constant multiplier (1/a?). 


(21) and, for the next eigenvalue ?, Rayleigh’s ratio px? is used. Values of the 
beam functions for different normal modes are obtained from tables of Character- 
istics Functions compiled by Young and Felgar’, and values of definite integrals 
are given in the Appendix. 


In Table II the lower and upper bounds are given for the first few mode shapes 
of skew plates for the three different sets of boundary conditions. The mean value 
is shown, and the percentage difference of the upper estimate from the mean. 
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Skew _ | 
angle | 
0° 1/2 1-328 
2/2 1:774 
15° 1/1 | 2:694 
7406 
30° 1/1 | 6-799 
45° 1/1 | 11-360 
(b) 
| | 0-772, 
| | | 1:125 
1-649 
| | 
30°C} | 5-490 
45° | 9-200 
1-783 
1/ | 1-057 
| 
15° 5-254 
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From the estimate of the lower limits it is apparent that with square plates very 
close bounds for the first few eigenvalues can be obtained by using only one term in 
the series representation of plate deflection. The “ brackets” in this case are very 
narrow and the results so obtained should for all practical purposes be equal to an 
exact solution. With skew plates such a simple approximation is not quite adequate 
and it becomes necessary to use more terms in the plate deflection series with 
increase in the order of the mode and the skew angle. For plates satisfying boundary 
conditions w=0 on the edges and having skew angles less than 30°, closer brackets 
to the natural frequencies can be obtained by using two or three terms in the series. 
But for angles greater than 30°, the estimates obtained are very poor; the widening 
of the inequality brackets is due to the fact that for higher skew angles the beam 
functions do not form a rapidly converging set. It is therefore necessary that, for 
large skew angles, either many more terms are used in the series representation of 
plate deflection, in which case the numerical computation becomes very laborious, 
or, if orthodox analysis permits, a more suitable function may be sought. 
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Appendix 
TABLE OF DEFINITE INTEGRALS 
The definite integrals which enter in the calculations are listed here. 


(i) Clamped-clamped beam 


a 
0 


8 


= (k,4- kin’) { k,a cot (k,a/2) k,a cot (k,,a/2) } (m, n even) 


8 kn? ky? {k a tan (k,,a/2) — k,a tan (k ( dd 
a? (k,* — m ma/4)— K,ala m,no ) 


= { 2 k,,a cot (k,,a/2)— {k,,a cot (k,,a/2)}? } : (m even) 


a~ 
| 
a. { {kna tan (k,,a/2)} 2k,,a tan (k,,a/2) \ i (m odd) 
where k,,a is the m positive root of the transcendental equation tan (k,,a/2) 


=(- 1)" tanh (k,,a/2) or of its equivalent cosh k,,a=sec k,,a. 
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(ii) Clamped-supported beam 


0 9 


{k atan k,a-—k,atank } ( ) 


= { k,,a tan k,,a — (k,,a)? } (m=n) 


where k,,a is the m‘ positive root of the transcendental equation tan k,,a=tanh k,,a. 
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The Behaviour of Damped Linear Systems in 
Steady Oscillation 


R. E. D. BISHOP 
ADDENDUM TO THE PAPER PUBLISHED IN MAY 1956 


In this paper (published in The Aeronautical Quarterly, Vol. VII, Part 2, p. 156, 
May 1956), a set of equations is derived giving the responses at the generalised co- 
ordinates g,, q., . . - Qn, Which are caused by a harmonic generalised force ®,e‘ 
corresponding to g,. They are equations (26), namely 


B B 
=i tz | 


§,-©? 


B ] 
| Be. 


This note is concerned with a sentence which follows these equations which says (of 
them) that “... there is no longer any assurance that the ratios between the constants 
mB,, in any column of partial fractions are the same for all co-ordinates g, at which 
excitation is applied, since the previous limiting condition cannot now be applied. Thus 
the shapes of the ‘ modes’ which are associated with the various columns depend upon 
the nature of the applied harmonic excitation.” 


The “ previous limiting condition’ referred to the conservative system and meant 
letting w approach a natural frequency at the same time as the amplitude of the 
excitation is diminished. In fact, this assurance can be obtained even for the damped 
system by letting 


—_ Ems ?, —> 0, 


where m=1, 2,....n. But, in using this limiting process, the reservation must be 
made that it has only mathematical (and not its previous physical) meaning. Thus it 
is not necessary to draw the conclusiofi as to the variation of the resonant modes with 
the point of excitation. 
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Suppose that the limiting process is used. A set of ratios is found, namely 


where the /’s are the amplitudes of the g’s. Now this ratio is that which is found to 
correspond with the root w? = €,, of the equation 


A=0, 


where J is the n'" order determinant formed from the coefficients of equations (17) of 
the paper. It follows that it is immaterial which co-ordinate qg, is used in establishing 
the ratios by the limiting process since there is only one set to be found from the 
equation 


This argument must be regarded as a purely mathematical one and, to see this, 
it is necessary to examine equations (17) of the paper. If the constants ®, are all zero 
in those equations, the relation A=0 is the condition that finite complex solutions for 
the W’s may exist. Now if the system is conservative, so that all the d,, terms are zero, 
then the equations (17) are those of free vibration and the solutions referred to determine 
the principal modes; one set of ratios exists for each root of A=0. But if the d,, are 
not zero, the homogeneous equations have no meaning, since the motion will not be 
harmonic and the d,, have meaning only when it is. 


It is commonly found that a modal shape of an aircraft which is found by a ground 
resonance test is dependent to some extent upon the point of excitation. For instance, 
the aircraft may be excited at a point on a wing at a resonant frequency and the 
corresponding modal shape determined by means of suitably-placed accelerometers. 
If, now, the excitation is applied athwartships with the same frequency at a point on the 
rear fuselage, the aircraft will again resonate. But discrepancies will be found between 
the second modal shape and the first which are such that the motion of the tail unit 
will be more pronounced in the former. 


This effect is due to the existence of damping forces. It can be explained in a 
rough-and-ready fashion by consideration of the problem as one of elastic wave 
propagation within the structure. It will also be seen that it is not inconsistent with 
the foregoing argument that the modal shape associated with the numerators of any 
column of fractions is the same no matter what is the driving point of the system 
concerned. 


When a particular column of fractions is brought to resonance, so that its contri- 
bution to the response predominates, other columns will account for some of the motion. 
The extent to which these extraneous modes modify the predominating mode depends 
upon (a) the damping forces present, (b) the proximity of other resonant frequencies 
and (c) the point of excitation. For, while the complex ratios between the numerators 
of other columns do not depend upon the point of excitation, the actual values of 
the numerators do. 
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